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1 Indroduction 

This review describes the processes of generation, evolution, and degradation 
of structures - linear, quasilinear, and nonlinear - and relaxed halos of dark 
matter (DM) in the Universe, based on theoretical developments and compar- 
isons with observations and numerical simulations. Special attention is paid 
to: (1) the characteristics of dark energy (DE) and methods of measuring 
it, (2) DE, which gravitationally affects the growth of cosmological inhomo- 
geneities and drastically influences the dynamics of DM structure formation, 
and (3) questions of normalizing density perturbation spectra on the basis of 
observations of the Universe's large-scale structure. We separately consider 
the current status of DM equilibrium halos (their internal structure, density 
distribution, rotation curves, etc.) in review [1]. 

Turning to the most challenging questions of physical cosmology, we do 
not attempt to highlight all aspects related to the formation and decay of 
large-scale structure in the Standard Model (in particular, we do not elabo- 
rate on the baryon history of the structure). Both reviews rely on the original 
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research of the authors and follow respective chapters of monograph [2]. Here, 
we propose a more complete account of the aforementioned topics, invoking 
new observations and comparisons to the theory, and avoiding detailed an- 
alytical manipulations whenever possible (all necessary theoretical results 
are formulated in the Appendices; the interested reader will find proofs in 
book 12]). 

2 What the structure is 

The Hubble flow (recession, rushing outward, expansion) of matter observed 
on large scales bears no structure: it is compatible with spatially homoge- 
neous and isotropic field of distributions of density, pressure, velocities, and 
other characteristics of matter. The structure is, by definition, inhomoge- 
neous and is associated with distortions of the Hubble flow that have evolved 
from initial seeding inhomogeneities of the flow metric under the action of 
gravitational field gradients, which influence the motion and distribution of 
matter in space (see Appendices A and B). Spatial gradients grow if scales 
decrease. Relatedly, one distinguishes between linear, quasilinear, and non- 
linear structures. 

By the large-scale structure of the Universe is meant the observed inhomo- 
geneous matter distribution which deterministically evolved from the initial 
small geometric scalar perturbations (the S-mode of cosmological perturba- 
tions) 'imprinted' in the gravitational potential of the Hubble flowQ The 
modern Universe exhibits a well-developed non-linear structure in the form 
of halos of galaxies, groups, and clusters at small scales {R < 10 Mpc), and 
shows a more regular, quasilinear distribution of matter on larger scales up 
to hundreds of megaparsecs, exemplified by superclusters and cosmological 
voids. There are numerous transition forms between nonlinear and quasilin- 
ear structures. 

Gravitationally bound halos are distributed nonuniformly in the Universe. 
They are frequently observed in sheet-like formations - 'walls' whose trans- 

^ These perturbations are also called adiabatic or growing adiabatic density perturba- 
tions. Theoretically, one may also conceive of primary perturbations coming from inho- 
mogeneities in the matter composition under invariable gravitational potential (so-called 
isometric perturbations). There are, however, no observational data (within the arrow- 
bars of the measurement) which indicate that a part of the initial conditions might be 
described by these perturbations. 
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verse size does not exceed 10 Mpc. A wall does not expand transversely and 
appears to be a nonlinear formation. It can be quasilinear along the other 
two more extended directions, continuing to expand along them. The walls 
themselves are inhomogeneous, and appear as a collection of elongated 'fil- 
aments' which may intersect, forming 'nodes' - the rich clusters of galaxies. 
The mean inhomogeneity scale in the Universe measures 10 Mpc (the density 
contrast variance in a sphere of this radius is close to 1, aio — 1), varying 
locally by increasing in the vicinity of clusters, and decreasing away from 
them. 

The observed structure is associated with the development of gravita- 
tional instability in nonrelativistic coUisionless matter. Initial quasiisotropic 
expansion of matter is accompanied by the development of anisotropy. In re- 
gions of augmented density, the expansion of gravitating matter slows down, 
comes to rest, and is superseded by a collapse. The initial stage of collapse 
proceeds mainly along one of three directions and leads to self-crossings (caus- 
tics) and the formation of one-dimensional oppositely directed flows. Later 
on, the regions of matter self-confined by its gravitational field relax, grad- 
ually acquiring a spherical shape and forming multistream systems trapped 
by gravity - the halos of DM. These processes, both well understood and 
studied, are corroborated by numerical experiments simulating billions of 
coUisionless gravitating particles (the A^-body simulations), and agree largely 
with observations (see, for example, Refs [3111]). 

Quasilinear structures (the density contrast Sr < 1 [see Eqn (62) in Sec- 
tion 13]) are subject to a rigorous analytical treatment since there is a small 
parameter enabling the machinery of perturbation theory in this case. There 
are two analytical approximations describing nonlinear formation - that of 
Zel'dovich [5] (exact in the one- dimensional case), and Press-Schechter [6] 
(exact in the spherical case). As indicated by numerical simulations, the 
Zel'dovich approximation describes well the large-scale matter distribution 
in regions where the collapse has just begun. In contrast, the Press-Schechter 
formalism pertains to fairly small scales and describes the distribution of viri- 
alized halos of DM. 
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3 Why galaxies form 



We observe the state of the Universe bilhons of years after the Big Bang. 
On large scales, the expansion of matter follows Bubble's law, which does 
not make a distinction between spatial points of the medium because the 
relative recessional velocity 5V of any two neighboring elements of matter is 
proportional to the proper distance 6r between them: 

5y = ^ = H5r . (1) 

The proportionality coefficient H = H{t) does not depend on the spatial 
coordinates or mutual location of these elementsU but depends on the proper 
physical time t. This law of matter expansion is preserved as a relict from 
the early history of the Universe, being, in essence, synonym of the notion of 
the Universe. It conceals all the information on the formation of the Universe 
and the seeds in its structure. 

Even though H is independent of spatial coordinates x, the metric de- 
scribed by Eqn ([1]) is homogeneous and isotropic only /oca//?/jfl i.e., it still 
depends on x = (i = 1,2,3). The Bubble flow ([1]) can be conceived of 
geometrically as a three-dimensional spatial (initially curved) hypersurface 
that uniformly stretches with time, preserving local isotropy at all spatial 
locations, with the factor of local stretching a{t) - ajj(x) , where a = a{t) and 
(x) are some smooth continuous functions of class C ^ . 

Indeed let us assume by definition 

Hi = aik a^^ , 

where aij and a^^ are direct and inverse symmetric positive definite matrices 
(the dot over a letter implies a partial derivative with respect to t). From 

^The modern value is Hq ~ 70 kms^^Mpc^^. 

■^Summing up vectors of distances and velocities on the hypersurface t = const, one 
readily concludes that if law ([1]) holds true for a particular observer, this same law is 
then also valid for all other points with just the same coefhcient H. This, however, does 
not prove the global homogeneity of space: we are dealing with a tautology, because an 
isotropic and homogeneous hypersurface is assumed a priori in this conjecture in the form 
of linear superposition of distances and velocities. Notice that the dependence of H on 
X can manifest itself on scales in excess of the size of the observable Universe. The real 
scale of the Friedmann world - the region of homogeneity where the linear superposition of 
velocities is valid - exceeds the Hubble radius [see Eqn (7) below]. Finding it is a question 
of the accuracy of observational data. 
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Eqn ([T]) we understand that they admit the factorization 



if/ = HSj : ttij = a ■ aij , a = exp (^J ^ dt^ , (2) 

and correspond to the interval squared between medium points: 

ds"^ = dt'^ — Sv'^ = dt^ — a^Qij dx^dx-' , (3) 

6t = 6ri = a ■ aij dx^ , gtj = gij{-^) = aikaijS''^ 

[the Kronecker symbol here follows from Eqn ([2])]. The converse statement 
is also true: differentiating 6r from Eqn ([3]) with respect to time we recover 
Eqn ([1]) with H = d/a. Relatedly, Eqns ([1]) and ([3]) are equivalent. They 
describe a homogeneous distribution of matter with a laminar Hubble flow 
such that the linear law of velocity superposition holds in the vicinity of any 
of its points. Passing from one spatial domain to another one, we have to 
redefine the rules of distance addition with regard to functions gij . 

The modern value of the scale factor is taken to be equal to unity: a = 
(1 + z)"^, where z is the redshift. In order to determine the boundaries 
of the Friedmann world, we expand gij in the Taylor series in some finite 
neighborhood of an arbitrary point: 

gij{x.) = Cij + ko CijkX^ + Cijkix'^x^ + . . . , (4) 

where the coefficients ko and Qj... depend on the selected point, which we have 
placed at X = in this case, and k^^ is the scale of variation of trace ■ju or 
the convergence radius fco|x| < 1 within which the first term Cij of the series 
exceeds the remaining sum. Since we are dealing with the form gijdx'^dx^ , 
we can always reduce Cij to the unity form Cij = 6ij by an appropriate choice 
of coordinates and leave only nonvanishing coefficients Cijk,,,- Expression (4) 
reduces then to the following: 

^?,,(x) = e-2''o(5,,.-26,,), (5) 

where the irreducible scalar go = and small tensor &q = ©jj(x) depend 
on X and are close to zero at x ~ (the subscript '0' on these functions 
implies the absence of the functional dependence on time) El 



*To be fully rigorous, we have to add to the right-hand side of Eqn (5) some scalar 
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According to Eqn (5), the quantity g(x) is defined up to an additive 
constant, and the difference in its values between distant points can be ar- 
bitrarily large. It is usually assumed that its mean value go in the observed 
domain equals zero. In that case, the deviations of g(x) from zero grow with 
the distance from the observer, and the upper bound on the Friedmann world 
is set by the condition 

\qiK')\r^l. (6) 

This size is a fortiori larger that the radius of external curvature {ko < Hq), 
since |g(x)| <^ 1 within the scale of observable cosmologyl^ 

qn = \q{H,')\ ^ 10-\ (7) 

Functions (5) encode information on the S and T modes of cosmological 
structure, which define the anisotropy of cosmic microwave background radi- 
ation (go and &o) and the seeds of galaxies (only go). It should be borne in 
mind that there are small corrections ~ 'Vq^/H to the metric (3) and that 
in this order the Hubble fiows are weakly distorted. The gradients of go may 

B^ij and vector terms depending on x (the comma in subscripts denotes a partial 
derivative over x). These terms, however, lack physical sense and can be removed through 
coordinate transformations. Vector modes, as well as decaying branches of scalar (S) and 
tensor (T) perturbation modes, are incompatible with isotropic expansion (2). In contrast, 
the growing branches of S and T modes of geometric inhomogeneities are preserved in the 
form (3) on large scales and do not violate Hubble's expansion law (2). Actually, these 
'frozen' gravitational potentials qo and ©o constitute the seeds of cosmological structure. 
The decaying modes had already faded out before the formation of the Hubble flow began 
and certainly had been negligibly small by the beginning of the galaxy formation epoch. 
^To derive relationship (7), one needs to make use of the Sachs- Wolfe formula [7] 

Stsw = qHz/5 ~ 10"^ 

and estimate qh with the help of spectral integral 

where qiiz = qkr^Ha is the large-scale Harrison-Zeldovich spectrum (8|f9|, and fccq = Ocq — 
0.01 Mpc^^ is the Hubble scale at the instant of time when radiation and matter densities 
are equal. Notice that relationship (6) is valid under the assumption of | ©0(^0^^) I < 1- 
The function ©0 carries information on cosmological gravitational waves and partly on 
the pre-inflational geometry of spacetime (the other part is hidden in the function go)- It 
is not, however, related to galaxies and, according to observations, is much smaller than 
(Jo- We do not consider it in detail for these reasons. 
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evolve and disrupt the laminar flow at small scales, leading to its breakup 
into self-gravitating nonlinear clumps of matter. 

We see that the early Universe is deterministic and all its motions can be 
considered in the framework of the Cauchy problem. By solving dynamical 
equations, we uncover the cause- consequence chain of events fully determined 
by the initial cosmological conditions (functions go and &o), which in fact 
gives rise to the cosmological time arrowj^ 

In summary, there are two regimes of matter organization, which are 
evolutionarily connected with each other and describe opposite processes: 
Hubble flows (large scales), and structure (small scales). From the GRT 
equations it follows that the metric (3)-(5) represents the leading term in the 
expansion of the exact solution in the small parameter 

f<l. (8) 

where /3 is the mean speed of sound in the medium, and k and H = aH = a 
are the spatial and Hubble frequencies, respectively. Of principal importance 
is the answer to the question of whether corrections to expression (3) grow 
or decay with time. The evolution leads to the breakup of existing Hubble 
flows in the first case, and to their creation in the latter. 

The answer to this question depends on the sign of function 7, where 

H dlnS 1 / 3p\ 
7 = -7F^ = --7T— = 1 + - ' (9) 



iJ2 dlna 2 V e 



7--^ = l + 7=2(l + - 



Indeed, in the first order in gradients of go we have [2] 

da 



g = go + y^p— ' (10) 

Vpec = -z^Vgo , b = gAqo , $ = 0go , (11) 

where g = g(t, x) is the curvature potential [TOl [H], is the gravitational 
potential, Vpec is the peculiar velocity of matter motion relative to Hubble 



^The time arrow can be violated in certain regions of spacetime where the relativistic 
effects are important (for example, in black holes or wormholes). 
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flow (1), and 5p and 5^=5 are the comoving perturbations of total pressure 
(p) and matter energy density (e) (see Appendices A and B). The growth 
factors V = z/(a), g = g{a), and (p = (f){a) depend only on time or the scale 
factor a{t) respectively: 



A critical regime for the growth functions is the linear cosmological ex- 
pansion for which 



In this case, u, g, and (j) are constant in time, and the evolution regimes 
mentioned above are separated: the Hubble flows are preserved at locations 
where they have existed, and no generation of new structures takes place. 

For a decelerated expansion (7 > 0), the functions z/(a) and g{a) grow 
with time, a gravitational instability takes place, the initially laminar medium 
flows become disturbed, and the conditions arise for inhomogeneous struc- 
ture formation on the side of small wavelengths (the function H monotonically 
decreases) . 

For an accelerated expansion (7 < 0), the functions u and g decay with 
time, and a new structure is not created but, in contrast, the build-up of 
Hubble flow and function g(x) over an increasing range of scales ko < k < H 
continues (the function H monotonically grows). 

One may conclude that gravity in equal degree spawns two dynamical 
properties: repulsion (the generation of Hubble flows), and attraction (the 
generation of structure). Which of them will prevail depends on the equation 
of matter state [see Eqn (9)]: the inflation (repulsion, the generation of the 
Hubble flow) is realized for {e + 3p) < 0, and the deceleration (attraction 
and the development of collapse) in the opposite case for {e + 3p) > 0. Being 
inherent to GRT, both inflation and collapse occur for rather general initial 
distributions and properties of matter and lead to the emergence of ordered 
geometrical configurations on various scales - the Hubble flows, and nonlinear 
halos of matter. 

The convergence of integral (10) at the lower limit assumes the dominance 
of initial adiabatic perturbations. If the isometric pressure scalar is neglected, 
the relationship between 6p and 6 takes a simple form: 6p = and the 
dynamics of scalar q obey the independent harmonic oscillator equation pUj 




(12) 



7 = 



a oct . 
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[TT] (see Appendix B) 



,,(,H.2^),-,^^^0, a^^^_^ (13) 

where G is the Newtonian constant of gravitation. For small velocities of 
sound 1/3 1 < 1, the leading solution q = go(x) continues into the causally con- 
nected domain [cf. expression (8)]. Spectral amplitudes of scalar curvature 
qo{k) in the post-recombination epoch and those of initial perturbations of 
the S-mode qk are linked linearly: 

go = T{k) ■ qk , (14) 

where T{k) is the transfer function of linear density perturbations that ac- 
counts for their evolution in the pre-galactic medium (for details, see Ref. [1]). 

4 Boundaries of homogeneity 

In the course of inflation in the early Universe, the stretching of the space-like 
hypersurface with the size of kQ^, already born by the preceding evolution, 
and the building-up of function g(x) from the side of small scales k ^ ko, 
which enter the zone k < H{t) from the microscopic region if H > took 
place. In this case, the growth of scales oca(t) takes over the Hubble radius 
H~^{t) (Fig. 1), all folds and irregularities of the initial hypersurface within 
k^^ are smoothed out, and the newly emerging small scales appear embedded 
into the already existing large-scale framework of the hypersurface being 
isotropizied. We can say that in the course of accelerated expansion of matter 
the minimum comoving scale of the Hubble flow decreases. 

The maximum size of Friedmann hypersurface k^^ is determined by the 
conditions at the beginning of inflation and, as follows from observations, 
exceeds the modern horizon (/cq < 2 ■ 10^^ Mpc~^, see Eqns (6) and (7)). 
The minimum size is linked to the end of the inflationary period of 
the Big Bang; it is certainly less than the size of the observable structure 
and compares to the wavelength of the background radiation quantum. In 
the process of post-inflationary decelerated matter expansion, the opposite 
process took place: the Hubble flow broke up from the side of small scales 
and its minimum size grew with time. 
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Figure 1: The Hubble radius (the hne with bends) and sizes of per- 
turbations oc a (inchned straight hues: 1 - a galaxy, 2 - a supercluster) as 
functions of a, from the Big Bang (B) to periods of radiation (r), DM(m), 
and DE (E). 

In order to determine the minimum scale of current Hubble flow (1), we 
write down the first terms of metric expansion in the Lagrangian coordinates 
(tc, x) (the label 'c' of comoving time is dropped where possible): 

ds^ = {1 - 26p) dt^ - 6r\ (15) 

Sr = 6ri = ae~''^{6ij - b go,ii) dx^ , (16) 

and rewrite Hubble's law in a more precise form [cf. Eqns (92) and (93) in 
Appendix A] 

where 

/dci 
h — . (18) 
a 

Neglecting the effective speed of sound, we have at the lower boundary ki of 
the Hubble flow: 

H 

6^1, Vpec ~ ^ ' ?0 ~ < 1 • (19) 

11 
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Figure 2: The spatially mean number of galaxies N{R) in a sphere of radius 
R (according to Ref. |12j). The accuracy of power-law asymptotics is about 
10%: 1 - oc R^, the uniform distribution of matter, and 2 - N oc R?, the 
nonlinear structure of the Universe. 

Recalling that the observed mean inhomogeneity scale ki ~ 0.1 Mpc~^ 
(Fig. 2), we obtain for the spectral amplitude of the curvature potential at 
the lower boundary of the current flow: 



in agreement with the required value of the transfer function T{ki) ~ 0.1 . 

Of interest is the question of determining experimentally the size k^^ of 
the Friedmann 'background' of our Universe. Doing so is possible, in prin- 
ciple. In fact, the scale k^^ need not be very large, because it is connected 
with the last period of accelerated expansion in the chain of inflationary 
stages of the Big Bangj^ and this period could be relatively short. In that 
case, k^^ can, in principle, exceed the radius of external curvature H^^ only 
slightly. In this situation, the total curvature potential qo{x.) generated to- 
ward the end of inflationary explosion and determining the geometry of the 
observable world would be composed of two components: a non- Gaussian, 

^Inflationary stages could have alternated with those dominated by matter which 
formed through the decay of intermediate short-lived inflatons. 



qoiki) -5-10 
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(20) 
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strongly-correlated smooth part of size kQ , linked to global pre-inflationary 
geometry, and small-scale (with respect to kQ^) ripples of a Gaussian field 
of inhomogeneities born in a quantum-gravitational way. The non-Gaussian 
component of density perturbations grows with an increase in the distance 
from the observer and manifests itself in the large-scale anisotropy of cosmic 
microwave back- ground. The Gaussian perturbations do not decay as the 
scale is reduced, and are responsible for the formation of galaxiesjf] 

The post-inflationary matter-dominated stage of decelerated expansion, 
accompanied by a reduction in the speed of sound, leads to the conditions for 
collapse at small scales, where a 'window' of gravitational instability opens. 
In contrast to the inflation, where initial conditions are forgotten, certain 
seed perturbations of curvature are required for the onset of collapse. They 
define domains of matter inflow and outflow. 

Thus, the scale interval /ci < <C fcm in the modern epoch belongs to 
nonlinear cosmological structures, and the undisrupted quasi-Hubble flow of 
matter, albeit already distorted, still persists at scales k < ki (in the mean 
in the Universe). 



5 The quasi- Friedmann equations 

A weakly inhomogeneous Universe is described by the generalized Friedmann 
equation [2] (see Appendix A) 

^v-^^-^, (21) 
in which the geometrical scalar variables 

i/, = i/.(t,x)^^ = i<^ 

(the semicolon indicates the covariant derivative) and b = b{t, x) correspond 
to the local Hubble and scale factors of medium volume expansion, e and 



^The resultant amplitude of quadrupole anisotropy of relic radiation can be both lower 
and higher than the value expected for this extrapolation of the short-wave spectrum. This 
effect could explain the low quadrupole anisotropy of the cosmic microwave background 
radiation (if, of course, this anomaly is confirmed by future observations and rigorous data 
analysis). 
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are the total energy density and matter four- velocity, the small function 

2 

x=x(t,x) = -Ag (22) 
3 

is the internal space curvature, and the dot above the symbols implies partial 
derivative with respect to the comoving time tc- The time evolution of matter 
density obeys the conservation law 

e + (e + p,) = , (23) 

where is the volume pressure [see Eqns (77) and (85)]. 

Equations (21) and (23) describing the evolution of the quasi-Friedmann 
Universe in geometrically invariant variables are valid for small spatial curva- 
ture |x| <C 1 and include the zero the first orders of magnitude in deviations 
of weakly inhomogeneous geometry from Friedmann's one. These equations 
have a scalar form, although the geometry described by them is anisotropic 
[see Eqn (89)]. 

In order to solve Eqn (21) in the volume factor b, the defining scalar of 
curvature q must be known. Assuming that the small function q is known to 
us, we seek the solution in the form 

H^ = Hc + ^ , 24 

3aH ' ^ ^ 

_ A$ _ AS 

^ ~ + 4^ ' P^-P-+ ' 
where = X{tc) are the background functions of comoving time, and $ 
and 5* are arbitrary small functions of all coordinates. Substitution into Eqns 
(21) and (23) results in the correct coupling between scalars $ and q [cf. Eqn 
(86)]: 

H f 

$ = — a{-fq-S)dt, (25) 



with the arbitrary function S being not linked to the Einstein equations. 

It should be emphasized that quasi-Friedmann Eqns (21) and (23) imply 
no constraints on the medium physical properties. Similarly to the original 
Friedmann equations, they link the spacetime curvature with total energy 
density, its time derivative, and the pressure of matter. As concerns the 
physical state of the matter, it needs to be considered only in the derivation 
of the equation of motion for the curvature scalar q (see Refs [2l [11] and 
Appendix B). 
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6 Dynamical properties of the structure 



A physical reason for the emergence of nonhnear structure - the formation of 
galaxies from small primary curvature perturbations - is the gravitational in- 
stability of dark matter, most vigorously developing in the post-recombinant 
period of DM dominance]^ Since the initial pressure in nonrelativistic mat- 
ter is low (relative particle velocities are close to zero at any point in space), 
the cold medium freely moves in own gravitational field of quasi-Hubble flow 
and the initially small peculiar velocities and the contrast of matter density 
grow with time on all scales. Consider at a greater length the dynamics of 
the quasi-Friedmann model at the linear stage of developing scalar inhomo- 
geneities. 

Let us turn to a simple model of the late Universe, which accounts only 
for nonrelativistic matter 'm' (with density inversely proportional to the local 
volume, Pm oc b^^) and dark energy 'E' of constant density (pe = const). In 
this case, the mean speed of sound in the post-recombinant epoch is equal to 
zero [see Eqn (95), (3 = 6p = 0] and, consequently, the scalar q and spatial 
curvature x do not depend on time, and equation (21) simplifies. Multiplying 
it by {B/He)'^ and recasting in terms of dimensionless variables, we arrive at 



f(6)-x(x), (26) 




where 



fib) ^ ^ (p. + pe) = + 6^) > 1 , (27) 

3 X 
q = qo{x) = - H^q , x = x(x) = — 2 = A g , (28) 
2 

i^E = Hq^/VL^ ^ (5 Gpc)~^ is the Hubble constant of dark energy, Cm = 
^m/^E — 0.4 is the constant coefficient, one sixth of which is linked to 



^It should be recalled that in the CSM the initial perturbations in composition are 
absent, the background curvature equals zero, and the parameters of the energy density 
of components take the following values: DE (JIe ~ 0.7), nonbaryonic DM(i7M ~ 0.25), 
baryons (ilb ~ 0.05), and radiation (fir — 10^''). Notice that the 10% accuracy level 
of today's observations does not yet enable distinguishing the cosmological constant and 
evolving DE. This witnesses in favor of its slow evolution and allows considering general 
models of DE in the form of expansions in terms of the small parameter |w + 1| ^ 1, where 
the cosmological constant is the leading term in the series [see Eqn (44), w = Pe/pe]- 
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baryons and the remaining 5/6 to DM (in this approximation, both com- 
ponents move together). Obviously, the function f{b) attains a minimum 
f . ~ 1 at h^^ ~ 1 7 

An arbitrary small function of spatial coordinates i< describes the local 
normalized space curvature. We are interested in domains with the positive 
right-hand side of equation (26): 

x(x)<l. (29) 

In these domains, the matter density decreases monotonically with time [see 
Eqn (42) for details]. They comprise both superclusters (x > 0) and cosmo- 
logical voids (x < 0). 

The volume and background scale factors coincide at points x = (the 
expansion anisotropy can be large in this case): 

b^a{t) = -^, H = hJ-^. (30) 

1 + z a 

where f — f{a) is the growth factor of the Hubble velocity component 

V^^ = ///ex. (31) 
In a general case, in the linear order in x we obtain 

b = a(^l-^gkj , 6^ = gk, (32) 

H,^H{l-^hk^, h=^^j^, (33) 

where 5m = 5 is the density perturbation, g = g{a) and u = p{a) are, 
respectively, the growth factors of density perturbations and peculiar velocity 
(Fig. 3): 



7i/2 1 / fda\ _ 3^^ 3He [da 

= a = — \ a — H I — , 1/ = -Hei^ = — — / — 

Cm Cm V JoHj' 2 2a2 H 
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1 2 3 4 5 

a 

Figure 3: Growth factors of density g{a) and peculiar velocity z/(a) 
For 

3a 

a < 1 : ^ - Y ' 

a > 1 : g = ^^ax ^ 1-56 , 

a = 1 : ^ ~ 1.13 , 

Equations (26)- (34) describe quasi- Hubble anisotropic flows with the ef- 
fective Hubble function which depends on the observer's location. In the 
modern epoch, the function v shows a wide maximum, which is an indica- 
tion of the period of most intense structure generation. The position of its 
maximum corresponds to 2; ~ 0.2, the 90% level of its maximum value of 
z/jnax — 0.68 is reached at a ~ 0.5 and 1.4, and the 50% level is reached at 
z ~ 0.1 and 4 . 

The current epoch is therefore that of maximum peculiar velocities, which 
will persist for the cosmological time [TH]. The function u will decrease 
twofold when the age of the Universe reaches 35 billion years. Only then will 
it be possible to speculate about the beginning of the epoch of the fading out 
of peculiar velocities in all space domains where x < 1 . 

The function h{a) determining the deviation of the local Hubble factor 
from its background value is plotted in Fig. 4. The maximum /imax — 0.65 



3 
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z/ ~ 0.67 - 



17 



0.0 4 1 1 1 1 1 1 ' 1 1 

1 2 3 4 5 

a 

Figure 4: Function h{a) setting the distortion of the Hubble flow [T3] . 

is reached at 2; ~ 0.4 , while the interval of the values of h > /imax/2 is within 
the limits a G (0.1; 1.8), which corresponds to the range of the Universe's 
age, from 0.6 to 22 billion years. Figure 4 explicitly demonstrates that our 
Universe is at the stage where Bubble's expansion law is maximally distorted 
[h{a = 1) ~ 0.87 ■ /imax — 0.57], and that the recovery of Hubble flows in 
quasilinear domains of space will take about ten billion years. 

We can conclude that the formation of large-scale structure in the Uni- 
verse lasts for the period from 1 to 22 billion years after the Big Bang. The 
stage of suppression of the Bubble flow inhomogeneities, caused by the grav- 
itational influence of DE, has not begun yet, although dark energy has been 
dominating in the density of matter for 3.5 billion years already. The delay 
hinges on the fact that the dynamic influence of DE on the structure gener- 
ation in the Universe has just begun. Indeed, the characteristic time of this 
influence exceeds the current age of the Universe and amounts to H^^^ ^ 17 
billion years. 



18 



7 Anisotropic cold flows 



The field of peculiar velocities can conveniently be described by passing to 
the Eulerian coordinates in which the gravitational field is locally isotropic 
in a linear order in x at any space point (see Appendix B). The relationship 
between the Eulerian (y) and Lagrangian (x) coordinates we are interested 
in has the form 

y = X + ^ S , (35) 

where S = S(x) = — V g is the vector of displacement for an element of 
the medium relative to its unperturbed position, k = — div S = A g . The 
Lagrangian coordinate x does not vary here with time along the trajectory of 
the medium element and coincides with the Eulerian coordinate y, as t — )■ 0. 
The growth factor for the displacement, g = g{a), is simultaneously that for 
the density perturbation (32). The fact that the matter displacement relative 
to the laminar Hubble fiow is factorized as a product of two functions (one 
dependent on time, and the other dependent on space coordinates) indicates 
that the rate of perturbation growth is the same for all wavelengths. 

The net displacement of medium elements with respect to the unper- 
turbed Hubble positions grows monotonically with time and amounts today 
to 14 Mpc. The mean net displacement will approach its limit value of about 
22 Mpc (see Fig. 3) in the future, provided that the DE density does not 
change. 

Using equation (35), we obtain the following representations for the in- 
terval on the Eulerian and Lagrangian grids: 

ds^ = (1 + 2$) dr^ - a^dy^ = dt^ - a^{5ij - 2^ g,^) dx'dx^ , (36) 

where t = t — avH-^q, and a = a{t) -(1 — g) = a(r) -(1 — $) is the local 
scale factor [2]. The function b{t,x) [see Eqns (21) and (32)] is proportional 
to the trace of the spatial part of the Lagrangian metric tensor, and the 
gravitational potential of density perturbations is equal to $ = 0.6 ■ (pq, 
where 

For a <^ 1, we get 0=1. The temporal factor = 0(a) of the gravitational 
potential decay under the action of DE is plotted in Fig. 5. The magnitude of 
can serve as a measure of DE dynamic infiuence on the structure generation. 
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Figure 5: Decay of the gravitational potential 0(a) of the large-scale structure 
of the Universe under the action of dark energy (see Ref. |2]). 

By the current epoch, the potential 0(a = 1) ~ 0.77 has already dropped by 
23% off its constant value in the matter-dominated phase. 

Interval (36) determines the physical Eulerian coordinate of medium el- 
ement r = ay ~ ay (cf. Eqn (105)). Differentiating r with respect to 
the proper time, we arrive at the following formula for the peculiar matter 
velocity: 



V 



pec 



- Hr = g'S = jjHeS. (38) 

Expression (38) coincides with the definition of 3-velocity as the spatial com- 
ponent of the 4- velocity of matter in the Eulerian reference frame: 

Vpec = TT- = -yH^yq = — . 

aay H 

Thus, the quantity u featuring in Eqn (33) is indeed the growth factor for 
the peculiar velocity. 

According to Eqn (38), the total velocity of matter is defined as 

V = Hr + Vpec • 

The first component Hr describes the Hubble velocity of the nonuniform 
Universe. It can also be split into two parts. One of them, V// = -ff x [see 
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Eqn (31)], is connected with the homogeneous component of density, and 
the other, gSS, is connected with the perturbed one. The growth factor 
for the second component is proportional to the growth factor g for density 
perturbations, whereas the pecuhar velocity is oc g'. 

We turn now to local flows of matter. In regions of the inhomogeneous 
Universe constrained by the condition (29), the flows of matter are described 
by the tensor field Hij = Hij{t,x) generalizing the function H{t) of the 
Friedmann model. Indeed, as follows from relationship (35), the coordinate 
distance between nearby points of the medium at time moment t is given by 

Syi = - gq,ij) Sx^ . (39) 

Differentiating the physical distance 5r = a5y over time, we obtain the field 
of pairwise velocities: 

6V, = ^{6r,) = H,^6r\ (40) 

Hij = HSij - gq,ij = H (Sij - hq,ij^ . 

The trace of the tensor field represents the volume Hubble function = 
Ha/ 3 [see Eqn (33)]; however, the tensor Hij itself is strongly anisotropic. 
The anisotropy of local expansion (variations of projections of Hij onto the 
radial directions emanating from the given point x) has the same order of 
magnitude as the deviations of H^ from the mean value of the Hubble pa- 
rameter H. At the boundary of quasilinear regions (29), these variations 
reach 100% (up to a point of stopped expansion in some directions). For 
example, in the vicinity of the Local Group, at a distance in excess of 
2 Mpc from its barycenter, the principal values of the Hubble tensor are 
Hij = diag(48, 62, 81 ) kmc'^Mpc"^ (Fig. 6). 

The field Hij describes regular cold flows of matter. It is noteworthy that 
Eqn (40) is valid under the assumption that the distance between galaxies 
is small, below the correlation radius of the two-point correlation function 
for the displacement vector. For different projections of this vector relative 
to the direction of Sy, the correlation radius varies from 15 to 50 Mpc. The 
deviations from the velocity field (40) grow with distance. The spectrum of 
cosmological velocity perturbations is shaped namely in this manner: it de- 
cays toward shorter wavelengths for k > 10~^ Mpc~^ (see, for detail, Ref. [1]). 
For this reason, random deviations from mean velocities (40) in this range 
grow with an increase in scale. 
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Figure 6: Radial velocities V{r) of galaxies in the vicinity of the Local Group 
according to data of Ref. [Hj. Different symbols correspond to different 
galaxy types, and distances to galaxies are given relative to the barycenter 
of the Local Group. 

It should be kept in mind that deviations from dependence (40) amount 
to ~ 40kmc~^, i.e., about 25% of the mean velocity (see Fig. 6) at a distance 
of ~ 2.5 Mpc from the barycenter of the Local Group. At the same time, the 
total peculiar velocity of the Local Group relative to the cosmic microwave 
background radiation comprises 600kmc~^. The scales of inhomogeneities 
responsible for so high velocities are in the interval from 15 to 70 Mpc. 

We see that the standard theory of the Universe's structure formation 
encounters no difficulties in explaining the observed motions of matter in 
quasilinear regions of the Universe (x < 1). The local flows are regular, 
smooth, and strongly correlated. The smallness of random galaxy velocity 
deviations from mean cold flow rates is explained by the shape of the initial 
spectrum of spatial density perturbations. These flows bear the quasi-Hubble 
character at small distances, preserving its main features: the flows are cold 
and radial, and the speed of galaxy recession is proportional to distance. 
However, the Hubble 'constant' depends on the observer position and the 
direction in space. We have already mentioned the neighborhood of the Local 
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Figure 7: Growth factors of Hubble /(a) (dashed line) and peculiar z/(a) 
(solid line) velocities of matter [2j; a = 1 at 2; = 0. 

Group as an example. The principal values of the Hubble tensor Hij at a 
distance of several megaparsecs from its barycenter are related as 3 : 4 : 5. 

As the radius grows, the deviations of galaxy velocities from the mean 
correlated flow rates increase too, beginning to saturate from distances of 
about 10 Mpc, which corresponds to the minimum correlation radius at which 
deviations of velocities reach the order of magnitude of the Hubble velocity 
proper. The deviations do not grow further, whereas the Hubble velocities 
continue to grow. At large distances, the flow of matter approaches the ideal 
Hubble law (l). 



8 How the collapsing flows form 

As shown earlier, dark energy, though not a component of structure, car- 
dinally influences the structure generation rate and the history of galaxy 
formation. Figure 7 illustrates the behavior of growth factors of Hubble and 
peculiar components of matter velocity at the quasilinear stage of evolution. 
Because of the peculiar velocity growth, the expansion in certain regions 
of the Universe turns into a collapse which starts in certain directions and 
proceeds then toward the formation of gravitationally confined halo systems 
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with various masses. The spectrum of initial perturbations is shaped so 
that increasingly larger masses may collapse at later time moments (see also 
Ref. P). 

If DE were absent, the Hubble velocity would continue to decrease with 
time, and the peculiar velocity would continue to grow. As a result, in 
any spatial domain where at least one of three principal values of tensor 
q^ij = diag(Ai, A2, A3) is positive [for example, Ai > 0, where Ai > A2 > A3; 
see Eqn (35)], sooner or later there will come an instant when the negative 
peculiar velocity in this direction becomes equal to Bubble's one and the 
expansion is halted [the principal value Hu turns to zero; see Eqn (40)]. In 
this case, the gravitational confinement of matter and its subsequent collapse 
ensue. However, because of the dynamical influence of DE not all domains 
succeed in evolving into the collapse stage: the Hubble velocity reaches a 
minimum at z ~ 0.7, and then increases sharply as oc a [see Eqn (31)], 
whereas the peculiar velocity saturates and then slowly decreases (see Fig. 
7). With account for DE influence, the process of galaxy formation resembles 
a fire burned out and lacking new firewood. 

Thus, the nonlinear structure has a chance to form in those spatial zones 
where condition (29) is violated, namely 

X = Ai + A2 + A3 > 1 . (41) 

Notwithstanding the threshold character of this inequality, it does not imply 
that regions occupied with nonlinear structures are topologically compact. 
Indeed, A3 (or A2,3) can appear to be negative, and then the expansion of 
matter will continue in that direction (or directions). 

How can the transition from the quasilinear stage to the collapse phase 
be described? The answer is surprisingly simple: solution (35) can analyt- 
ically be continued up to the first self-crossing {5yi = 0) and even further 
by 'matching' the arising multistream flows across caustics (the Zel'dovich 
approximation [5]). 

This analytical approximation leads to correct qualitative conclusions, 
while quantitative deviations from the real evolution (until the first self- 
crossing) do not exceed 20-30% according to different criteria. The success 
of the approach, which is so simple, hinges on the potential q being small 
under cosmological conditions (the formation of black holes is exponentially 
suppressed). Moreover, the solution (35) proves to be exact in the nonlinear 
one- dimensional case, and namely this variant of initial collapse development 
is most typical in the Universe. 
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Relatedly, when computing the quantitative characteristics of a newly 
forming nonhnear structure - distribution functions for cosmological voids, 
superclusters, filaments, etc., mean distances between the walls, nodes, etc., 
correlation functions, evolution based on redshifts, and so on - one can rely 
on linear perturbations and their Gaussian statistics because solution (35) is 
naturally set by the field of small initial density inhomogeneities. 

We give below two examples. First of all, we can exactly specify those 
domains where caustics fail to form, i.e., the expansion continues forever: 

Ai<^-aL-0.6. (42) 

As long as density perturbations remain small, the two conditions (29) and 
(42) are equivalent. Condition (29), however, by no means ensures the ab- 
sence of collapse: if Ai > 0.6, this direction, even if only in the future, will 
turn to the collapse. In a similar way, we can refine the density field. Prom 
relationship (35) we obtain the comoving matter density 

Pm(r/,x) = 



a^dei{5yi/5xj) a\l - g X,) {1 - g A2) (1 - ^ A3) 

Formula (32) follows obviously from the last one for ^ <^ 1. 

We will return to nonlinear structures in the subsequent sections, and 
now proceed with the measurements of DE. 



9 How to measure dark energy 

DE can be measured only by means of observational cosmology. Its detection 
in the laboratory (similarly, for example, to attempts at laboratory detection 
of dark matter) seems to be implausible because this ever-penetrating sub- 
stance (field) only weakly interacts with all objects, be they interiors of stars, 
compact objects, the early Universe, or others, during the entire world history 
known to us. Only indirect measurement methods are possible. Luckily, 
they have already proven their efficiency and led to impressive results - the 
discovery of DE proper. 

What does it mean to measure DE? The answer to this question is hinted 
at by cosmology. The known constraints on the main DE parameter, viz. 

\l + w\< 0.1 (43) 
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(see, for example, Ref. [13]) point to the constancy of function w = Pe/^e ~ 
const both in time and in space. Condition (43) simphfies the approach to DE 
detection to some extent. The DE clustering effects are of small significance 
and escape detection at the current level of technological development. 

Today, we can only discuss actual experiments dealing with the depen- 
dence of function w on time: w{a). Taking into account the slowness of 
evolution, this function can be replaced by the set of constant coefficients c„ 
of its Taylor expansion around —1 : 

w{a) = -1 + Co + CiAa + . . . + Cn^^^ + ■ ■ ■ , (44) 

where Aa = 1 — a = z{l + z)~'^. At the current level of knowledge, the 
physical information on DE available to us is stored in these coefficients. 
From inequality (43) we find | c„| < 0.1. Certainly, first of all, a pending 
question is the measurement of cq and the leading expansion coefficients. 
Theoretical predictions, based on connections of c„ with physical parameters 
of models {(po, fn, and others), in no way discriminate their range of values. 
Hence, the question hangs solely on the possibility of experimental assessment 
of the function w{a). 

We single out three main ways of measuring w{a): structural, dynamical, 
and geometrical. Let us consider briefly each of them. We have already dis- 
cussed the structural method earlier. Without being a part of the structure, 
DE cardinally influences the galaxy formation rate. It is the structural argu- 
ment that has led to the discovery of the DE phenomenon. The achievements 
of all of 20th-century astronomy established that only a small fraction of mass 
enters the large-scale structure of the Universe. The remaining largest part 
is, hence, contained in the form of unstructured DE. Figures 3-5 and 7 de- 
pict the growth curves for the seeds of structure in the case where all c„ = 0. 
For c„ different from zero, the evolution curves will deform, as can easily 
be computed with the help of the CRT. The point, therefore, lies in max- 
imally precise observational measurements of real growth functions and in 
determining the parameters c„ with their help. 

How can the growth curves be measured accurately? Unfortunately, one 
can hardly rely on the traditional astronomical methods of observing the com- 
pact sources at different wavelengths (stars, including supernovae, quasars, 
galaxies, clusters, and so forth) and assessing based on these data the quan- 
titative characteristics of inhomogeneous matter distributions in space and 
time. The obstacle here is uncontrolled nonlinear effects of coupling between 
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density distributions of matter and light. They, eventually, make impossible 
the recovery of mass distribution, based on the luminosity of its baryonic 
component, with necessary accuracy. 

The true breakthrough in the structural method of detection of DE prop- 
erties can be expected from dynamical measurements which are sensitive 
to gradients of total gravitational potential. An example can be furnished 
by statistical measurements of background galaxies weakly lensed by nearby 
structures, but even here barriers are met. The point is that the retrieved 
surface density of total mass (for a given redshift), which is the final product 
of the method, is still insufficient for a successful comparison with theory. No- 
tably, we are not in a position to compute the contribution from baryons with 
necessary accuracy because of their complex interaction with light (shock 
waves, stellar formation, supernovae, dissipation, ultraviolet background, 
ionization, cooling, and so on). 

In today's discussion of maximally precise methods of measuring DE, the 
case in point is largely the exploration of linear and quasilinear structures 
because, in this case, there is a well understood and fully controllable theory. 
An example of a breakthrough in this field is furnished by measurements 
of anisotropy and polarization of cosmic microwave background radiation 
(CMBR) that has led to the creation of CSM. A similar breakthrough is also 
possible by further developing this model based on accurate measurements 
of cosmological parameters, including those of DE. 

This breakthrough can be provided by any statistical measurements of 
large-scale peculiar velocities of matter, inhomogeneities of gravitational po- 
tential, and distributions of matter density. In the first case, we have to do 
with the Doppler effects pertaining to the motion of matter [161 [U] 5 and the 
measurement of the field of peculiar velocities based on proper motions of 
galaxies on the celestial sphere [I8j. In the second case, it is the retrieval of 
gravitational potential with the help of the weak lensing effect, but this time 
at larger scales, where self-crossing of streams of matter is absent and the 
contribution from baryons is easily assessed [19]. In the third case, we are 
dealing with quantitative distribution of structures as a function of redshifts, 
and the detection of baryonic acoustic modulation of the density perturbation 
spectrum (see, for example, Ref. |20j). 

Staying on this path, we arrive at the dynamical method of measuring DE 
by virtue of the Sachs- Wolfe integral effect. The influence of DE induces an 
additional source of the CMBR anisotropy, related to the decay of the linear 
gravitational fleld in the current epoch (see Fig. 5). By measuring this effect. 
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Figure 8: Superimposed portions of a cosmological microwave background 
map, which are projected onto 50 voids (a) and 50 superclusters (b) ~ 100 
Mpc (4°) in size, based on results of Ref. [21]. The colorbar of the anisotropy 
of cosmic microwave background radiation (the insert at the top) is in mi- 
crokelvins. 

one can exactly determine the properties of DE, as has already been done 
when determining the parameters of the early Universe from measurements 
of the anisotropy parameters of CMBR in the recombination epoch. Quanti- 
tatively, here we need to raise the measurement accuracy only by one order 
of magnitude. 

The dynamical integral effect predicts the existence of a cross-correlation 
between two celestial maps - those of CMBR anisotropy, and the large-scale 
structure of the Universe. Figure 8 shows superimposed regions of CMBR 
maps, projecting onto 50 voids (Fig. 8a) and 50 superclusters (Fig. 8b) 
set apart from the galaxy distribution in the SDSS catalogue. The mean 
size of voids/superclusters is about 100 Mpc. Summation of respective re- 
gions with reduced and increased densities of matter is performed to amplify 
the cross-correlation signal against the background of a random field of pri- 
mary perturbations. We observe an apparent valley in the temperature of 
relic radiation (a cold spot) in directions towards voids, and a hot spot of 
the appropriate angular size in the direction towards superclusters, in full 
agreement with theoretical predictions. 

These are only the first results which demonstrate the correspondence 
with what we already know well from the structural studies. Here, however. 
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the precise cosmological information is encoded; it can be deciphered in fu- 
ture observations and used to construct a more precise model of the Universe. 
Reahzation of this wonderful opportunity requires two exact maps of linear 
cosmological perturbations: the anisotropy of CMBR, and the gravitational 
potential of the large-scale structure. All this is possible at the already avail- 
able technological level, but hinges on building new-generation telescopes, 
both ground-based and operating in space. 

10 How to measure the Universe 

Impressive successes in exploring the Universe rely on differential and dy- 
namical measurements of gradients of gravitational potential, peculiar ve- 
locities, and matter density perturbations by methods of observational cos- 
mology. There are also, however, direct astronomical methods of measuring 
the zeroth-order geometry bypassing the structure. For brevity, we will call 
them geometrical. These methods are advantageous at large scales because 
the Hubble flows build up with distance, while deviations from them decrease. 

Geometrical tests deal with measurements of distances and times between 
events. In the zeroth order, these intervals are controlled by the function 
a{t) = (1 + 2;)^-'^ and its derivatives. Both the properties of DE and other 
parameters of matter can be reconstructed by the measured cosmological 
functions H{z) and 7 (-2). 

There are two classical astronomical methods of determining geometrical 
sizes: the measurement of radial distances or angular scales as functions of z. 
In the first case, we are dealing with the dependence 'apparent magnitude- 
redshift', m = m{z), and in the second, the 'angular size-redshift', 9 = 
0{z). Other geometrical tests exist, too (for example, counting the number 
of galaxies inside spheres of a given radius; see Fig. 2), but in this review we 
consider in more details only the two mentioned above. 

The relationship m = m{z) is called the Hubble diagram in astronomy. 
It allows determining the scale factor a{t) of the Universe if the absolute 
luminosity is known for the objects being analyzed. Indeed, (a~^— 1) is the 
redshift of the source, and {to — t) is the distance to it in light years measured 
along the light cone of the past. Astronomers frequently use the luminosity 
distance r, which defines the fiux J-' of radiation from a point source of 
luminosity Lq recorded on Earth, as if both the source and observer were in 
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the Euclidean space: 

Instead of J-" and Lq, one may use the apparent m and absolute M stellar 
magnitudes, which are related to the distance r between the source and the 
observer as 

m = M + 5 (Igr [pc] - 1) . (45) 

In order to estimate the absolute magnitude of the source luminosity, one 
needs to understand the source structure and the way it shine, i.e., to have 
a predictive theory. As a rule, this task is not solvable with the necessary 
accuracy, as we are dealing with a complex gasdynamic nonlinear system with 
many parameters, be it a galaxy or cluster of galaxies, a quasar, or a star. 
For this reason, one frequently resorts to empirical relationships that offer an 
estimate of object luminosity based on its other observable characteristics. 
These phenomenological relationships are built and validated on the nearest 
sources and then used to calibrate distant objects of the same class. It is 
implicitly assumed that near and distant sources of the same class are alike 
(the standard candle hypothesis). 

While different objects have been proposed to play the role of the standard 
candle over years0 this research cannot be used in rigorous cosmology as 
long as there is no theory for a system of that complexity, the theory based 
on numerical simulations and taking into account numerous, as yet little- 
studied, factors which pertain not only to the physical nature of objects (for 
example, star formation, nonstationarity, chemical composition, and nuclear 
fusion), but also to their environment, conditions of light propagation, and so 
on. For such data, a probability will always exist that we are dealing not with 
the effects of geometry or composition of the Universe, but rather with the 
internal properties of sources and their evolution (uncontrolled systematics). 

While a standard candle is required for making use of Hubble diagrams, 
one needs to know the physical size of the observed object (the standard 
rod hypothesis) in order to measure distances by its angular size. Let the 
physical size of an object be d, and its angular size 6. The distance based on 

^"By way of example, the brightest cD galaxies in centers of rich clusters were used 
in the 1970-1980s as the standard candle. Currently, supernovae, gamma bursts, X-ray 
clusters of galaxies and other bright sources that yield to observations at large distances 
pretend on this role. 
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the angular size, D, is then defined as 



Making use of the interval of the Priedmann geometry, we get a relationship 
between the distances introduced above: 

r = {1 + zf D = {1 + z)R, (46) 

where R is the geodesic comoving distance to the object: 

R , , - f dz 

R=Jf^-Vo-viz), R^H.j-. (47) 

(?7 is the conformal time). Apparently, the cosmological function R{z) de- 
pends only on the properties of matter, not on Hq. 

With the help of formula (46), any diagram m{z) can formally be recast 
as the dependence 0{z), and vice versa. This is in no way surprising, as 
we are exploring unified geometry, even if resorting to various observational 
tests. In summary, the recipe for measuring the Universe is simple: construct 
diagrams 6{z) and 0{z), or m{z) and rn{z), or others. If the physical size 
of the source or its absolute stellar magnitude are known, we will manage 
to determine the distance to it and, hence, the dimensions of the Universe. 
This will enable us to reconstruct the functions H{z), 7(2;), along with other 
functions of the scale factor, and, in turn, determine the composition and 
properties of matter and their evolution with time. 

The realization of this program requires highly accurate measurements 
and knowledge of the nature of the objects under study. In so doing, these 
objects can be very diverse. Figure 9 presents a cumulative diagram of 9{z) 
for several geometrical tests reduced to the standard size of d = 9.5 pc. 
Shown are the results of measuring distances based on the luminosity of 
distant la supernovae (points to the left of the minimum), the angular size 
of ultracompact radio sources (bases of jets from active galactic nuclei: six 
crosses in the lower part of the curve), and the anisotropy of CMBR (spec- 
trum acoustic oscillations: the cross in the upper- right part). The dashed 
fine corresponds to the CSM. 

The most reliable point of observational data in this diagram is the right 
cross determined by the angular size of the acoustic horizon at the moment of 
recombination ^^rec = 1100. The points and crosses found from observations of 
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Figure 9: Diagram of 6{z) based on three geometrical tests reduced to the 
standard size d = 9.5 pc (according to Ref. [22]). 

supernovae and ultracompact radio sources rely, respectively, on the standard 
candle and size hypotheses. There are certain physical grounds for using 
them. In the first case, we are dealing with the thermonuclear explosion of 
a binary white dwarf with total mass in excess of the Chandrasekhar limit 
setting the maximum luminosity of a supernova la. The second case pertains 
to the limiting mass of a central accreting black hole (~ 10^° M©, Mq is the 
mass of the Sun), which defines the 'standard size' of radio bright quasars. 
Note, however, that the available models of supernova explosions and active 
quasars do not meet for the moment the requirements placed on accurate 
geometrical measurements. Rather, the opposite approach can be beneficial: 
by knowing the answer (the parameters of CSM), we can inquire into the 
proper cosmological evolution and the nature of compact sources of that 
type. 

Figure 10 demonstrates the image of an ultracompact radio source with 
minimum resolution on the order of a parsec, attained with the help of a 
terrestrial radiointerferometric network. The base of a jet can be seen well, 
but the internal structure of the source itself is unresolved. This example 
illustrates that exploration of compact objects, which lays the basis for di- 
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Figure 10: Ultracompact central radio source in Cygnus A (adapted from 
Ref. |23]). Lyxs is the distance expressed in light years. 

rect geometrical measurements in the Universe, calls for the development of 
new techniques. Their realization is possible by making use of cosmic radio 
interferometry. 

With the help of space telescopes, the interferometric baseline can be 
increased and, respectively, its resolution can be improved. Such instruments 
not only are capable of measuring distant compact objects but will also 
detect time changes in their angular sizes and redshifts, accompanying the 
expansion of the Universe. The last two factors are directly linked to the 
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Figure 11: Trajectories of light rays in the Friedmann model in conformal 
coordinates (r/, x) from the source {x = R) to the observer (x = 0) (see 
Ref. [2]). 

Hubble function^ 

0_ H{z) z _ H{z) 

m-^°"rT^' ^^^^ 

This opens new prospects for determining Hq and H{z) at cosmological dis- 
tances. 

Notice that if the DE density were constant, the redshift of the sources of 
the Hubble flow located at zq ~ 2.3 would not change: z{z = zq) = 0. As this 
takes place, the redshift of more distant objects {z > zq) decreases with time, 
and increases [z{z < zq) > 0] for less distant objects. The boundary z = zq 
is rather sensitive to the evolution of DE. Its detection will help answering 
the question of how the DE density varied in the past. 

Determining the evolution effects (48) at cosmological distances will re- 
quire several years of monitoring at angular resolution on the order of mi- 



These relationships follow upon differentiating the functions 

^ d , , a(n) 

^ - and 1 + z - ^" 



Rai-q-R) a{r]-R) 

over 77 for constant d and R, where 77 is the observer time, and (77 — R) is the conformal 
time of the source (Fig. 11). 
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Figure 12: The plane {R,y) observer O and the jet node (the arrowhead) 
moving with velocity v relative to the source S pj. 

croarcseconds, which is technically challenging. However, even the first cos- 
mic interferometers will be in a position to measure distances to nearby galax- 
ies by their proper motion relative to the cosmic microwave background [T8] . 
Moreover, one may measure distances by observing superluminal motion of 
the nodes of jets in active galactic nuclei, directed at a small angle to the 
line of sight. Such measurements are already being conducted on Earth (see, 
e.g., Ref. [21]) but can be essentially improved with the help of cosmic inter- 
ferometers. 

Let us clarify this effect in Fig. 12. Let a be a small angle between 
the motion direction of a jet node (the arrowhead) and the line connecting 
the observer O with the source S (quasar), and rj and x = {R,y,0) be 
the conformal coordinates of the node [R is the geodesic distance from the 
observer to the jet; see Eqn (47)). The observation time for the jet node 
moving with the intrinsic velocity v is given by 

dr = drj — dR = (1 — v cos a) d-q , 

and the apparent velocity of its displacement in the observer tangent plane 
is defined as 

dy vsina 

vj_ = — = < vF . (49) 

dr 1 — V cos a 

This function passes through a maximum at cos a = v, equal to the jet 
gamma-factor F = (1 — v^) for the velocity v = |v| approaching the 
speed of light. 

The jet displacement on the celestial sphere takes then the form 

• v^ 0.3 / F \ _i 
fc'^ = — - < I — arcms year , (50) 



R - R V20, 

which constitutes milliarcseconds per year for the superluminal velocities 
of bursts from quasars observed in reality (Fig. 13). By having extensive 
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Figure 13: Superluminal ejections from quasars (see Ref. [25]): (a) the jet 
transverse velocity vj_(z), and (b) angular displacement 6±{z). 

statistics of such objects, we will be able to determine the geometrical char- 
acteristics R{z) (47) using the upper envelope of function (50). 

We offered several examples illustrating the potential of direct assessment 
of geometry and composition of the Universe from observations of nonlinear 
objects. They can be augmented by observations of bright galaxies, compact 
groups of galaxies, rich X-ray clusters, gamma bursts, binary quasars, and 
other active systems. Exploring the Universe with their assistance remains 
not as efficient as that using quasilinear systems and based on measurements 
of large-scale structures, microwave background, and nonrelativistic matter. 
We may hope that the development of theory and numerical experiment will 
make the geometrical methods competitive in accurately solving the tasks of 
observational cosmology. 
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11 Formation of galaxies in the inhomoge- 
neous Universe 



How do the gravitationally bound halos of matter, reached hydrostatic equi- 
hbrium, emerge out of the quasi- Hubble AovJB and how is the mass dis- 
tribution of these nonhnear objects hnked with the initial field of density 
perturbations? 

In contrast to the large-scale structure bound up with relatively small 
matter density variations and staying at the beginning of its nonlinear de- 
velopment, the galaxies have evolved from high density peaks typical for the 
field of small-scale perturbations (less than several megaparsecs) and have 
already passed through the period of nonlinear relaxation. The reason for 
that is the amplitude and shape of the initial density perturbation spectrum 
which increases toward short wavelengths (see Ref. [1]). This property of 
spectrum facilitates early fragmentation of matter into 'halo-blocks' of small 
mass and their subsequent agglomeration in more massive halos of galaxies, 
groups, and clusters. 

The hierarchical formation of nonlinear structures from small to large 
masses finds support in numerous observations and in numerical modeling. 
Important observational argument in favor of the sequential formation of 
halo systems is the absence of a large number of far emission Li/a lines. The 
explanation is that massive systems, in their past, did not undergo a single 
powerful burst of star formation involving all halo gas, but formed by way of 
subsequent merges of numerous small-mass blocks that passed through their 
star formation phases at different times. 

By virtue of the existing properties of the spatial S spectrum, the observ- 
able structure of the Universe possesses the following important features. 

— The formation of galaxies and groups is largely completed, they contain 
the dominant part of DM. The process of galactic cluster formation still 
continues. 

— Halos and their characteristics are distributed non-uniformly in space: 
their number and masses are modulated by the large-scale structure of the 
Universe. 

— The first feature offers physical motivation for establishing a simple 
linkage between emerging nonlinear halos and peaks in linear seed density 

^^We also call the equilibrium halos relaxed or virialized systems, with the understanding 
that it is the virialization over velocities, not energies. 
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perturbations, based on the model of a quasispherical collapse [6]. The second 
feature invites us to apply the Press-Schechter method to finite-size domains 
and, in this way, to couple the mass function of halos with the large-scale 
structure of the Universe. 

Before discussing theoretical predictions related to the halo mass func- 
tions and their validation against observations and numerical experiments, 
we remind the reader of the main characteristics and elements of the large- 
scale structure of the Universe. As we have already stressed, turning to large 
scales moves us to the initial stages of structure formation and the first phases 
of collapse, which are predominantly one-dimensional. Long-term observa- 
tions of the spatial distribution of gravitationally bound virialized DM halos 
on scales of up to 300 Mpc demonstrate that the large-scale structure of the 
Universe possesses the following typical elements: 

— filaments — linear structures of varying richnes^ with lengths of up 
to tens of megaparsecs, and thickness of several megaparsecs; 

— walls — flat formations reaching hundreds of megaparsecs over elon- 
gated directions and several megaparsecs in thickness, and filled with fila- 
ments; 

— nodes — rich clusters of galaxies severalmegaparsecs in size occurring 
at intersections of filaments; 

— superclusters — extensive domains of space from tens to hundreds of 
megaparsecs in size filled with walls, filaments, and nodes; 

— cosmological void^^ — extensive domains reaching several hundred 
megaparsecs, where bright galaxies and galaxy clusters are absent. 

The last two structural types are quasilinear because of their large size. 
Under the action of gravity, matter moves towards superclusters which ex- 
pand more slowly than the mean Hubble flow. The matter 'flows off' the 
voids in different directions, and their local volume grows faster than the 
mean. The mean displacement of matter elements for the entire history of 
the Universe is smaller than the sizes of these structures and makes up 14 
Mpc. The voids occupy the largest part of the Universe (more than 50%) and 

^■^The richness of an element of a large-scale structure is defined by the number of 
galaxies composing it. 

-'^'''The term cosmological 'void' reflects the fact that on celestial maps some regions 
seems to be empty. In reality, these cosmological empty domains can contain a significant 
(yet not dominant) fraction of the Universe's matter, i.e., they are not really empty. That 
is why, instead of referring to cosmological 'empty domains', the term 'voids' is frequently 
used. 
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their forms in contours of low density are close to spherical, whereas super- 
clusters are dense formations with a flattened shape. Both types of structures 
can naturally be characterized by the mean substance density within a given 
volume. It is smaller than the Universe's mean value in voids, and larger in 
superclusters. 

Observational assessment of the parameters of voids and superclusters is 
ambiguous. The reason for the blurring of voids and superclusters recon- 
structed from the spatial light distribution resides in a nontrivial physical 
coupling between luminous matter and the DM composing the backbone of 
the structure. The resolution of these formations through the threshold den- 
sity contrast of the matter contained in them compared to the mean level is 
more rigorous and accurate. We define voids (and superclusters) as three- 
dimensional spatial domains with a given negative (positive) level of matter 
density contrast. When the absolute value of threshold contrast decreases, 
the sizes of structures grow, and their number decreases. For p < pc, we 
obtain a single void filling the whole Universe. 

An important characteristic of voids and superclusters is the mass func- 
tion of the halos populating them. Observations indicate that the cosmolog- 
ical voids nearest to us are impoverished in normal and dwarf galaxies, but 
weakly emitting galaxies and gaseous clouds are encountered. They are ob- 
served both through the ultraviolet absorption lines in the spectra of quasars 
or bright galaxies placed behind them and the absorption line of neutral hy- 
drogen (21 cm) in the spectra of distant radio sources. Galaxies in voids 
are of reduced luminosity on average, compared to those in superclusters, 
though the slopes of their luminosity functions differ but slightly from each 
other |26] . 

The theory predicts that voids and superclusters contain weak galaxies 
of small mass that have not passed the stage of active stellar formation, as 
well as primordial objects of the hierarchical clustering model with masses 
up to 10 ^M© and less. Such objects are incapable of holding gas, so that 
stellar formation in them is practically impossible. With growth in mea- 
surement sensitivity, observers are succeeding in finding previously unknown 
weak galaxies with masses of 10^ — 10^ Mq in the Local Group and the near- 
est neighborhood [H]. Their number agrees well with the predictions of the 
GSM, yet the question of finding lessmassive galaxies remains open to date. 

The minimum cutoff scale for the spectrum of initial density perturbations 
is linked to the physics of DM and is still unknown. The region of small 
masses that yields to exploring by means of observational cosmology (for 
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example, by observing forest of absorption lines in the spectra of distant 
quasars) can be linked to the temperature T of the Universe at an instant of 
time when the comoving size i? of a given mass M/j coincides with Bubble's 
radius of the Universe: 

Numerical modeling of processes pertaining to the formation of structure 
from cold matter has led to substantial advances in exploring the mass func- 
tions of nonlinear halos in different domains of the Universe. Starting simu- 
lations of the formation process at 2; ~ 20 with seed masses M ~ 10^ Mq, to 
2; ~ one manages to obtain a spatial distribution of gravitationally bound 
systems with typical masses in the range ~ 10^^ — 10^^ Mq (see Refs [27] 
- which is similar to the distribution observed. In this case, all dark 

matter proves to be involved in the virialized halos of all sorts of masses, 
while the mean matter density inside the halos at the instant of their for- 
mation exceeds the mean density in the Universe by a factor of ~ 200. This 
corresponds to the predictions of quasispherical collapse. 



12 Mass function of relaxed halos 

The Press-Schechter approximation relies on two assumptions. 

(1) The linear perturbation field is Gaussian: the volume fraction of the 
Universe where the density contrast x), smoothed over a sphere of ra- 
dius R [see Eqns (61), (62)], exceeds a certain threshold value 6c is given 
by 

f{M,z) = ^ e-^"/'dv', (51) 

where 

v = v{M,z)^^ (52) 

is a monotonically increasing function of argument M = = Airp^R^/S, 
and 

aR{z) = g{z) ■ Or (53) 

is the variance of density contrast in a sphere of radius R (see Ref. [1] for 
more detail). The connection of the normalized factor of density perturbation 
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growth with those introduced earher is defined as 



'^'^ = 5(1) = 1(1) ■ '''' 

(2) The matter in density peaks collapses and stays confined in gravita- 
tionally bound objects: there is such a value 6c whereat the mass fraction of 
matter in the Universe!^ -^(> z) = 2/(M, z), turns out to be involved in 
virialized halos with individual masses in excess of M by time moment z: 

oo 

F{> M,z) =2f{M,z) = — [ ndM . (55) 

Pm J 

M 

where 

n = n{M,z) = ^fM (56) 

is the differential mass function (the mean number density of halos with 
masses M in the interval dM ~ M), and 



din (Xji 



f{u) (57) 



dlnM 

is the fraction of matter that resides in halos with masses M by time moment 
z. 

A comparison of expressions (51) and (55) allows us to determine the 
mean mass function of nonlinear halos in the Universe: 

i{u)=fps^^iye-'''/\ (58) 

The normalization condition of the function fM is written down as 

' dM 

„ f«ir = > (89) 



^^The numerical coefficient 2 implies that in the course of hierarchical clustering all the 
cold substance proves to be caught in virialized halos of all possible masses: F{> 0, z) = 1 
for any z. The argument M is more convenient than R because M is preserved both 
for linear perturbations and for halo objects. If the component of massive neutrinos is 
included in the matter, we need to take into account the dependence of growth factor (53) 
on the spatial scale. 
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and it implies that for any z all the dark matter is confined in nonlinear 
halos.We can only inquire about the distribution of these halos over masses 
and about the typical halo objects for a given z, containing a dominant part 
of the Universe's matter. For large 2, all the matter is decomposed into small 
masses, but the share of large halos grows with time. The mass of typical 
halos, M^,, at any epoch is defined by the condition 

v{M,,z) = 1. 

For the modern Universe, one has M^, ~ 10^^ Mq. For small M, the function 
iu ~ ^ 1. This implies that halos of small mass are seized by more 
massive ones in the course of hierarchical clustering, i.e., they leave the set 
of M objects and are taken into account only in newly forming halos of larger 
mass. 

The approximation of linear density contrast in the CSM at the instant of 
halo formation in the framework of homogeneous coUaps^ gives the thresh- 
old contrast value 5c = 1-675, which only weakly depends on variations of 
cosmological parameters within their error bounds |33] . 

Long-term investigations have shown that the results obtained with the 
help of this method agree well with the numerical iV-body experiment. Nev- 
ertheless, the basic formalism does not include mechanisms of merging and 
tidal breakup of objects. Hence, as a rule, the Press-Schechter method is ap- 
plied to analyzing the spatial distribution of objects with a sufficiently high 
mass (such as clusters of galaxies), for which these effects are insignificant. 

In order to apply this method in precise cosmology, account must be taken 
of corrections for the nonsphericity of collapse. Analytical estimates show 
how the function i{y) is modified in this case. For constructing an accurate 
analytical approximation, we need to use a more general empirical formula 
enabling a nonspherical correction (the Sheth-Tormen approximation [3^): 

i{y) = fsT ^ fo (1 + e-'"/' , (60) 

^^In a simple cycloid model of a collapsing dust sphere a = uq (1 — cos{y/>cri)), the time 
moment z of halo formation coincides with that of collapse, t{z) = j^adri = 27rao/v^- 
Taking into account that the radius of a virialized halo is half that of the radius enclosing 
the same mass at the moment the cycloid stops, we conclude that the density of the halo 
at the instant of its formation exceeds that of a collapsing sphere at the instant of its 
stopping by a factor of 8. Hence, inter alia, a useful estimate follows for the halo radius 
at the halo formation instant - it is approximately six times smaller than the radius of a 
sphere enclosing the mass M in the unperturbed Universe. Speaking figuratively, a halo 
at its inception is a fragment of the unperturbed Universe compressed sixfold in scale. 
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Figure 14: Functions f(z/) in the Press-Schechter (PS, dashed hne) and Sheth- 
Tormen (ST, sohd hne) approximations of halo mass functions (see Ref. [2]). 



where u = v/vq., and find numericaUy the correction factors 
fo = 0.32, z/o = l-2, p = 0.6 

(fo is the normahzation coefficient obtained from Eqn (59)). 

Such a corrected approximation, accounting for the Gaussian perturba- 
tion field and the collapse of density peaks, ideally describes the results of 
numerical simulations and is widely used in modern cosmology (Fig. 14). 
Apparently, the Press-Schechter approximation somewhat over estimates the 
number of gravitationally bound halos with masses near the characteris- 
tic M^, {v = 1), but underestimates the number of massive halos with 
M > 5 X 10^^ Mq. The intersection of two functions happens a.t u = 2.3, 
which corresponds to the mass of order ~ 10^^ Mq in the current epoch (the 
precise value depends on the normalization of the spectrum; see Section 14). 

In its standard form, the Press-Schechter formalism offers expressions 
only for the mean mass function of virialized halos, and misses such large- 
scale inhomogeneities as voids or superclusters, in which the local density pi 
of matter on the scale L differs from the background one, pm- In the next 
section we will obtain mass functions for halos populating the inhomogeneous 
Universe \35\ . 
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13 Modulation of galaxies by the large-scale 
structure 



Local extended domains in the Universe can conveniently be characterized 
by the mean matter density inside a sphere of comoving radius L centered 
at a point r: 

PL = PL(r, z) = j p{r', z) WL{\r - r'|) dv' (61) 

where 

WL{r) 



47rL3 [ , r>L. 

Let us introduce the parameters of mean density contrast and local matter 
density in the domain L: 

5l = ^-i, nL = n^^. (62) 

Pm Pm 

Adding up expression (62) with Qe, we obtain the mean density of matter 
in the domain L. Subtracting from it the equality Qra + = 1 for a flat 
Universe, we obtain the parameter A^, characterizing the mean curvature of 
domain L: 

AL = nL-n^ = n^5L . (63) 

Distinct from R, the free parameter L is assumed to exceed the inhomo- 
geneity scale {L > 10 Mpc), so that spatial domains with < p^ (voids 'v', 
Fig. 15) and Pl > Pm (superclusters 's') are still in the phase of quasi-Hubble 
expansion {\Sl\ < 1). As L — )■ oo, p^ tends to the mean matter density in 
the Universe, pm = Pmi^)- 

In order to construct the local halo mass function we shoTild: 

— resolve the field of linear perturbations of the density contrast 6r into 
a large-scale background 6l and a small-scale part which characterizes 
the seed perturbations leading to halo formation in a given spatial domain 
of size L > R (see Fig. 15): 

Sr = Sl + Sr\l , (64) 

— describe a quasispherical collapse on the background pi instead of the 
homogeneous background of density pm- 
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Figure 15: Spatial field of density contrast ^^(r) (solid line) as the sum of 
long-wave background 5l (dashed line) and short-wave residual 5^1^ (dotted 
line) [35j; 5c is the threshold contrast. 

From a physical viewpoint, we are considering the influence of the local 
large-scale structure of the Universe on the local process of hierarchical clus- 
tering. Assuming the scale of the structure is large (L > _R), its influence 
on the collapse can be represented in the form of a power series in spherical 
harmonics. The monopole term in this series is the main factor influencing 
the local statistics of virialized halos, which lays the basis of this method. 
The dynamics of collapse depend only on the mass of the substance confined 
within a sphere of radius i?; however, the regions where the density contrast 
exceeds the threshold one are less abundant in voids than in superclusters. 
Namely this information is contained in the background offsetting the Gaus- 
sian function 5r\l. At this stage, we ignore the dipole (density gradient), 
quadrupole (tidal forces), and other terms in the structure expansion, which 
are linked to the motion of the halo as a whole and the nonsphericity of col- 
lapse. Their influence can be accounted for as numerical corrections to the 
halo mass function [see Eqn (60)]. 

The dispersion 5r\l of a linear Gaussian field describing the local proper- 
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ties of density perturbations for R < L takes the fornix 

oo 

0-^2^ = j P{k)\W{kR) -W{kL)\^k'^dk. (65) 



Apparently, aR\L — )• for M = Mr — )■ (the mass in the sphere of radius 
L). 

Carrying out computations similar to those in Section 12, we obtain the 
halo mass function M G (0 , Ml) in the domain L |35] : 

nL = nL{M)^^^M\L. (66) 

fjWIL = 

where the function f(z/) coincides with Eqn (58) (or Eqn (60) in the case of 
correction for nonsphericity) , and Vl follows from Eqn (52), with replaced 
by (Tr|l, and 5^ by 

The function = i>i{M,z) grows monotonically with the mass growth 
and diverges when M — )■ Ml. The threshold density contrast 6c\l depends 
on the local density of background. When computing 6c\l it should be borne 
in mind that the domains of superclusters with 6l > (or voids with 6l < 0) 
expand more slowly (faster) than the homogeneous background. A smaller 
(larger) threshold contrast 6c\l with respect to the actual background will as- 
sure the collapse of a given domain of increased density up to a time moment 
z: 

Sc\L = Sc-6l- (67) 

This formula is valid in the linear order in 5^. Numerical correction is needed 
to account for nonlinear contributions. 

Figure 16 presents halo mass functions MuL/pm at z = in a void, 
supercluster, and flat domain of space with dimensions L = 140 Mpc and 
different values of mean density contrast 6l- As can be seen, small-mass 
halos form the main population of zones with a reduced matter density. 

For example, the fractions of substance in a void with 6l = —1/2) existing 
in the form of primary isolated blocks - halos with masses 10'' Mq, 10^ Mq, 

^^Here, the field crosscorrelation between Sf^L and 6l is neglected. It is significant only 
at i? ~ L. 
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Figure 16: The DM fraction (1 + 6l) iM\L in the halo of masses M (in the 
interval SM = M) in the domains of the modern Universe with the size of 
L = 140 Mpc and density contrast 6l (see Ref. [35]). The curves correspond 
to a supercluster (dashed, 6l = 1/2), flat background (solid, 6l = 0), and 
void (dotted, 5^ = -1/2). 



and 10^^ Mq - constitute about 1, 10, and 30%, respectively. It is easier, 
therefore, to discover low-mass primary objects of DM, which are not confined 
in more massive halo systems, in voids than in other parts of the Universe. 

The solid line for spatially flat domains of the CSM in Fig. 16 character- 
izes the modern Universe as a whole. We see that the range 5-10^^— 5-10^'^ Mq 
of halo masses (massive galaxies and groups of galaxies) hosts 30% of the to- 
tal DM in the Universe, and massive halos with M > 5 x 10^^ M© (rich 
clusters of galaxies) include no more than 10% of its nonrelativistic matter 

(fM < 0.1). 

Thus, the process of hierarchical gravitational clustering, collecting the 
dominant part of low-mass objects into more massive systems with M < 
5 ■ 10^^ Mq and observed as a phase of the active merging of galaxies, is 
completed in the Universe on the whole. For more massive formations with 
M > 10^'^ Mq this process, seen as the interaction of groups and clusters, 
still continues and is far from completion!^ 



^The patchy structure of most galaxy clusters in X-ray spectral region offers an addi- 



47 



If we consider separate domains of space, the coalescence of galaxies and 
interaction of galaxies there followed their own ways in the past and continue 
to do so. Indeed, clustering of halos in voids fully ceases with time [see the 
limit mass Ml in Eqns (65) and (66)]. In contrast, in 'fiat' zones and in 
regions of augmented density, the clustering and coalescence are progressing 
and gradually shift toward larger masses. We may argue that the activity 
in the form of coalescences, accompanying the process of galaxy formation 
everywhere in the Universe since its beginning, degenerates with time, con- 
tinuing in flattened superclusters and breaking down in voids. The process 
resembles a burning out fire. 

As a result, the morphology of galaxies proves to be different at various 
locations in the Universe. Notably, finding isolated galaxies which lack ad- 
jacent partners of comparable mass is, for instance, more probable in voids 
and close to their boundaries than in superclusters. In contrast, interacting 
galaxies are more common in superclusters, beginning with pairs and ending 
with large associations or clusters. 

Summarizing, we can conjecture that the theory considering the forma- 
tion of the nonlinear structure from DM agrees well with observations and 
numerical experiments and has a predictive skill valuable for further research. 
There still remain some problems related to the evolution of the baryon com- 
ponent. Solving them hinges on further development of computer facilities, 
observational techniques, and data processing methods. 

The halo mass functions in voids and superclusters differ substantially. 
In particular, the spatial density of massive (M ~ lO^^M©) halos in voids 
is several times less than in superclusters, and the distinction in mass func- 
tions is even larger for larger masses. In regions of galaxy clustering, the 
observational statistics worsen: 



where M* ~ 4 • 10^*^ Mq on the average in the Universe. For this reason, the 
comparison of theory against experiment relies on integral functions of halo 
masses, possessing more extensive data statistics. 

The integral mass function describes the full number of halos Nl in the 

tional argument in favor of this statement. 




(68) 
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Figure 17: Integral mass functions Nl{> M) of equilibrium halos in the 
range 10^^— lO^^M© for the same domains of the Universe as in Fig. 16. 

domain L (reduced to a unit volume) with masses in excess of a given one: 

Ml 

Nl{>M) = j ULiM)^. (69) 

M 

The differences between the integral functions in voids and superclusters 
diverge exponentially with increasing M, because the characteristic values 
M£ of halo masses are different there (Fig. 17). In the region of small 
masses (dwarf galaxies), the integral functions are hardly distinguishable, 
and the f^iL functions appear more appropriate (see Fig. 16). The actual 
distinction between mass populations of voids and superclusters is stronger 
because of the nonlinear galaxy coalescence effects and tidal destruction of 
galaxies, which partly suppress the small-mass part of the mass function in 
superclusters. 

The main population of voids is primordial objects (POs) with small 
masses. The share of the matter in voids in the form of isolated galaxies 
with M < 10^ Mq reaches 10%, whereas it does not exceed 1-2% in other 
domains of the Universe. Globular clusters in the central regions of massive 
galaxies can be stellar remnants of POs: they could have formed in molecu- 
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lar clouds of the first small-mass halos and then accumulated in the central 
regions of massive galaxies in the process of hierarchical clustering. 

The primordial objects in voids can manifest themselves as weak dwarf 
galaxies or in the form of Lya absorption systems. Indeed, we can estimate 
their spatial density and dimension, assuming that their formation epoch 
Zq ~ 10 and mean mass M ~ 10''' Mq, as 

npo ~ ^ ~ lO^Mpc"^, rpo ~ ^n^Q^ {1 + zq)~^ ^ Ikpc. 

As can be seen, these primary formations are rather loose, so that lensing 
effects on them are unimportant. The possibility of discovering them de- 
pends on their subsequent evolution. One may hope that such objects have 
retained their gaseous component in voids and that their column hydrogen 
density appears sufficient for the emergence of weak Lya absorption lines. 
The probability of picking such objects up in the line of sight in a void of 
size L ~ 100 Mpc is on the order of unity: n^or^QL ~ 1. 

14 Normalization of scalar perturbations 

Data on the anisotropy of the CMBR give exact information on the ampli- 
tude and shape of the initial S spectrum in the early Universe and on certain 
parameters of composition (densities of nonrelativistic components, and the 
curvature of space). However, to trace the transformation of initial inhomo- 
geneities into the observed galaxies, we also need to know the growth factor 
for density perturbations in the late Universe, which depends on the Hubble 
radius and DE (see Sections 6 and 7). The number and mass distribution of 
gravitationally bound relaxed halos, formed through the development of the 
gravitational instability of DM, depend exponentially on the amplitude of 
initial scalar inhomogeneities in the curvature and growth factor of density 
perturbations (see Sections 12 and 13). Owing to this dependence, we have 
at our disposal a sensitive test enabling us to determine the amplitude and 
shape of the S spectrum together with the most important supplementary pa- 
rameters of the Universe's composition from the quantitative characteristics 
of the structure. 

This test is frequently referred to as the normalization of the spectrum of 
cosmological density perturbations. For a given observed mass distribution 
of gravitationally bound DM halos, the amplitude of the power spectrum 
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is a function of cosmological parameters (fim? ^e? and others) within their 
variabihty bounds admitted by the data accuracy. 

For historical reasons, the normahzation 'to sigma 8' - the dispersion 
of density contrast within a sphere of radius 8h~^ (in this section, h = 
ifo/100kms~^Mpc~^), which represents the integral function of the density 
perturbation spectrum - is used most widely. We denote this dispersion as 
cTii here to emphasize its actual value: 8h'^ ^ 11 Mpc. The sphere encloses 
the mass of an unperturbed Universe, Mn ~ 2-10^^ Mq, which is close to that 
of a typical cluster of galaxies. This normalization test is therefore ideally 
oriented to using observational data on the abundance of galactic clusters. 
Moreover, the theoretical analysis of halo mass functions in this mass range 
relies on an elaborate analytical formalism (see Sections 12, 13). 

The magnitude of dispersion an, which ensures the observed spatial num- 
ber density of clusters, essentially depends on the total matter density Q^a in 
the Universe, whereas changes in all other cosmological parameters affect an 
only within 10-20% of its magnitude. We will demonstrate below the tech- 
nique of normalizing the density perturbation spectrum to the abundance of 
galaxy clusters at z = 0, which satisfies the current requirements of precise 
cosmology [36j. 

All cosmological parameters can be spread between two levels according 
to their impact on the quantity an: 

— (the first level); 

— Qe, h, n, fj,, Qb, We, and others (the second level), where = Q^/fl^ 
is the fraction of matter in the form of massive neutrinos. Among the free pa- 
rameters of the second level, the first four are statistically significant (within 
intervals of their variability). The rest are fixed for simplicity {Qi^h"^ = 0.023, 
and We = —I)- 

A class of models to be utilized represents a rather advanced variant of 
the CSM extension at the modern knowledge level: it includes the nonzero 
spatial curvature {fl^ = 1 — f2m — ^^e), hot, cold, and baryonic components 
of matter {Q^ = + Qu + ^b), the cosmological constant, and nonflat (but 
power-law with the slope n) spectra of the S mode. In this class of models, 
the exact dependence of the threshold density contrast 6c on and Qe in 
the curved Universe is also taken into account [33] . 

The subsequent study is split into two stages. At the first stage, an 
optimum value of an with accompanying individual error was computed 
through the comparison of theoretical and observational differential functions 
of cluster masses for every realization of the extended model (with its own 
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Figure 18: The dispersion aui^flj^) of density contrast in a sphere of 11 
Mpc in radius for the CSM The vertical lines show the accuracy of 

observational optical data on the abundance of near clusters of galaxies (a 
confidence level of 95%). 

set of cosmo logical parameters) @ 

At the second stage, all the computed values of an (they are tens of 
thousands in number, and each has its own error) were fitted by the approx- 
imating dependence of the form ^6] 

an = n^'+^' '^-+^^ [A + (70) 

+ As {n^ - A(i){l- Ajh- Asn- Ag Q] , 

The optimal values of parameters Ai{i = 1,2, ...,9) were determined by 
the Levenberg-Marquardt method of minimization. Notice the important 
distinction of this stage from the previous one. In the first case, n(M) was a 
function of a single variable (mass M). In the second case, the function an 
depended on multiple variables {Qui, ^e, h, n, and f^). 

""^^The approach presented here can work with arbitrary observational data pertaining 
to virial masses of galaxy clusters. In the particular case considered here the results of 
optical observations of the speeds of galaxies in the nearest 150 clusters with a median 
value of redshift z ~ 0.05 (see Refs [37l[38]) were utilized. 
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Figure 19: Projection of the space of cosmological parameters in the extended 
CSM onto the plane ((Tii,r2m) for the Sheth-Tormen (bars labeled with 1) and 
Press-Schechter (labeled with 2) approximations (from Ref. |36]). 

Figure 18 illustrates the second phase of normalization. It shows the 
slice of the cosmological parameter space of the extended model in the two- 
dimensional plane of main parameters (an, fim). In this case, the standard 
values were taken for other cosmological parameters (^e = 1 — = h = 0.7 , 
n = 1, f,y = 0). The plot confirms the expected, strong dependence of ctu on 

The experience of working with different data indicates that the value 
obtained for an depends on the selection of observations. We have already 
mentioned that this is, unfortunately, a common story accompanying work 
with nonlinear objects (clusters of galaxies in this case). This should be 
borne in mind when mentioning the 'standard values' of, first of all, such 
parameters as an, Hq, and f^m (or Qe)- They, to a larger degree than 
any other quantities, depend on observations of astronomical objects (stars, 
quasars, galaxies, and clusters). There is only one way out: improvement in 
the quality of observational data and control of systematic effects related to 
the evolution of the baryonic matter component. 

We now turn to other parameters. The fact that their contribution to an 
is statistically significant is vividly illustrated, besides by formulas [see Eqn 
(71) below], by Fig. 19. All the space of models is projected there onto the 
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plane (cth, flm)- As can be seen, the height of bars an for various values 
of parameters h, n, and f^, is comparable to the size of statistical error, 
and exceeds it for some sets of cosmological parameters. This emphasizes 
the importance of a thorough account of the contributions from all model 
parameters mentioned above. 

An interesting result, also seen from Fig. 19, is the increasing influence 
of the second-level parameters as the values of fl^a reduce. Since the interest 
of modern cosmology is focused namely on this domain, this once again 
confirms the necessity of simultaneous account of all cosmological parameters 
when solving cosmological tasks at a high level of accuracy. The scatter of 
(Til upon variations of second-level parameters looks tighter for the Sheth- 
Tormen approximation, which says something about its advantage: it ensures 
the most accurate power spectrum normalization in work with mass functions 
of galaxy clusters (see, for instance, Ref. [39]). 

We specify the form of formula (70) for the normalization of the scalar per- 
turbation spectrum as obtained in Ref. |36] , which compared the abundance 
of optical galaxy clusters with the theoretical halo mass function derived in 
the Sheth-Tormen approximation in the extended CSM: 

^^^Q 0.54 + 0.15 Qni- 0.04 He _ (^^1) 

-0.2 {n,n- 0.75) (1 - 0.2 /i - 0.2 n + 0.8 /^) = 0.53 ± 0.08. 

Here, the left-hand side contains combinations of cosmological parame- 
ters, whereas all errors are mapped into the right-hand side and attributed 
to the normalization (a confidence level of 95%). 

This example illustrates the specifics of cosmological research with respect 
to observational data: 

— the need for multidimensional spaces of parameters to determine the 
main model characteristics, 

— the need for numerous multivariate observational tests to lift the prob- 
lem of parameter degeneration. 

As follows from the example above, formula (71) alone is insufficient to 
normalize the S spectrum amplitude itself: we only know the relationship 
between the parameters, and need additional data to solve the normalization 
task. Note also an interesting feature of this formula: it has the form of a 
'plane' for the second-level parameters. The coefficients of this plane may be 
refined by more complete data of future observations. 
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15 Conclusions 



Summing up, we can argue that the standard theory for the structure for- 
mation in the Universe in the framework of the CSM does not encounter 
principal problems in explaining the observed matter motions and distribu- 
tions. The theory of nonlinear structure formation hinges on the dominance 
of DM. It agrees well with observations and numerical experiments and has 
predictive skills for further research. There remain complex problems per- 
taining to the nonlinear dynamics of DM and the evolution of the baryonic 
component (we consider some of them in Ref. [1]). Solving the problems 
listed above and those related to them will benefit from the development of 
computer facilities, new methods of observation and data processing, and the 
refinement of cosmological parameters aimed at further improvement of the 
theory and extension of the CSM. 
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16 Appendices 

A. Quasi-Hubble flows in the GRT 

Appendices A, B, and C present the main features of inhomogeneous cos- 
mological models in which the Friedmann symmetry is fully broken, yet the 
departure from it remains small [2]. Such models, deviating only slightly 
from the spatially fiat Friedmann model, are called here weakly inhomoge- 
neous or quasi- Friedmann, the flows spawning them, quasi-Hubble, and the 
geometrical departures proper, considered in the linear order of smallness, 
cosmological perturbations. 

One of the most important geometrical characteristics in the GRT is the 
Einstein tensor, or the stress-energy tensor of matter T^. Energy and mo- 
mentum are carried in spacetime along timelike world lines - flow trajectories. 
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A vector tangent to them, u'^, is the eigenvector of tensor T^: 

T;u^^ = eu\ (72) 

and its eigenvalue e = pc^ is the scalar of total comoving density@ The unit 
vector {u^u'^ = 1) describes the transport velocity of total flow energy. 
By definition, one has 

where = 5^ — u^^u'^ is the projective tensor orthogonal to u^. Accordingly, 
the symmetric tensor PavT^ describes the pressure structure. 

In the Friedmann Universe, solution (72) is unique for all forms of mat- 
ter: the components of a medium move along the vector (in Friedmann's 
coordinates = 6q). The stress-energy tensor of matter has a universal 
form which contains the only 'new' scalar p = p{t), the effective pressure in 
the medium: 

T|^ = {^ + P)u^^u'' -p6l^ = eUf,u'' -ppI^. (73) 

For weak perturbations of the Friedmann group, different medium com- 
ponents move in different directions, but deviations of their velocities from 
the mean velocity ■u'^ are small. Under this condition, solution (72) is unique 
and we obtain the general form of in the quasi-Friedmann geometry: 

T; = ie + p) u,u^ -p6;-s;: = e u,u^ - {pp; + s;) , (74) 

which contains, in addition to the form (73), a small stress tensor orthog- 
onal to the flow {u^s^ = 0)- As we see, the total tensor of quasi-Hubble 
flow pressure comprises two terms: the diagonal tensor p, and the anisotropic 
stress tensor s^. The total pressure can be determined by the trace of the 
pressure tensor: 

Ptot = ^P^%.=P+^<. (75) 

In a weakly inhomogeneous model, as in the Friedmann model, the rules 
of linear superposition are valid for matter components interacting only grav- 
itationally: 

e = S , p = J:pm, = S/„,M^ , (76) 

^''in the energy units of measurement adopted by us, energy and mass densities coincide 
(c — 1). Their different notations e and p, respectively, are only important for dimension- 
ality recovery and passages to the limit. For example, Gp does not contain the physical 
constant ~ speed of light in vacuum - and possesses the same dimensionality as H^. 
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where scalars of partial contributions of matter components have the form 

fm^'-^^^^, ^fm = l, m = l,2,...,M. 
e+p 

For general nonlinear interactions, the resolution of densities and pressures 
into components is ambiguous; however, notions of partial enthalpies = 
£m + Pm and component velocities u^, as well as total e and p, are preserved. 
Notice that these simple linear superpositions of medium component veloci- 
ties are valid for terms of zeroth and first orders of smallness with respect to 
deviations of velocities u!^ from the mean one u^. 

Projecting the Bianchi identities T^.^ = onto and orthogonal di- 
rections, we arrive at equations for the energy and momentum of the total 
flow: 

u''e,, + 3H^{e+ptot) = 0, (77) 

a^ = u''u^.^^u'' {u^^^-u^^^) = ^^^^ , (78) 

e + p 

where = u^.^^/'i is the Hubble factor of medium local volume expansion, 
and is the flow acceleration {ai,u^ = 0). To make the resolution of pressure 
into isotropic and anisotropic parts unique, we demand that the scalar part 
of s'^^ be of zero divergence [see Eqn (81)]. 

In the class of quasi- Friedmann coordinate splits (t,x), where the three- 
dimensional flow velocity is always small {D and v are small functions): 

u^^il + D,Vi), (79) 

the decomposition of geometrical variables into the background and pertur- 
bation parts yields 

g^^u^gil^ + h^u, (80) 

= diag (0 , -a%,) , lh,.-(^ ^2(^4^/| B,.^ ) , 

£ = + Se, p = p^^^ + 5p , (81) 
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where sym denotes a symmetric matrix, and the letter F (discarded below 
when possible) labels the Friedmann variables. They are functions of t and 
obey the equations 

= 47rG£(^) , + 3H{6^^^ + p^^^) = . (82) 

The scalar sector of perturbations is fully described by four gravitational (A, 
B, C, D) and four material potentials (f , fe, and 5* for arbitrary physical 
fields!). The scalar S of anisotropic pressure is gauge invariant, while the 
remaining seven functions are not Hll 

The gauge invariant variables for the dimensionless velocity potential and 
perturbations of density and pressure in the matter have the form 

q = A + Hv, 5 = ^, 5p = ^. (83) 

e+p' ^ e+p ^ ^ 

where the Lagrangian variables are 

6ec = Se — ev = e — Ec , (84) 

Spc = 5p-pv = p-pc = Ptot-Pv (85) 

The scalar of a volume pressure p^ = pc + s'^/3 and background functions 
Xc = X^{tc) describe distributions on spacelike hypersurfaces of constant 
comoving time tc = t + v. 

The field q = q{t,x.) plays a central role in the description of density 
perturbations [TOl H]. It has a double physical sense: on the one hand, it 



^""^It should be reminded that the vortex and tensor parts of sj^ are Unked with V and 
T perturbation modes and are not considered here. Decompositions (80) and (81) of 
geometrical objects into a background and perturbations are not unique: under a small 
coordinate transformation a;'' .t^ — we will obtain a new background (the same 
background functions, but for other time t) and new perturbations, but the total geometry 
will be preserved. Expanding a small arbitrary vector = (X, a^Y^i) in two potentials 
X and Y, we get the following gauge transformations for scalar variables: 

A- HX , B ^ B + Y , C^C + X + a^y, 
D^D + X, v^v + X, t^t-X. 
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is the dimensionless potential of the total matter velocity (laboratory inter- 
pretation), while, on the other hand, it is the potential of spatial curvature 
(cosmological interpretation). From the Einstein equations in the lowest or- 
der one derives the relationship between the potentials of the S mode and 
the field q: 

^ = 47rGfee, fbq-S)^, (86) 

a J H 

y-C + a'B^j^, v-D^j^, (87) 

= I , (88) 

where q = g — $ is the potential of the peculiar velocity of matter [see Eqn 
(101) below]. Hence, notably, having chosen v — C = 0, we obtain the 
metrics in the Lagrangian orthogonal reference frame {tc,^): 

ds^ = (1 - 2Sp) dtl - {5ij - 2B, ij) dx'dx^ , (89) 

where 

u^^il + Sp)S^, -^^/q^' a = a-(l-A) = ae-(l-g) 

is the scalar scale factor. The proper time ds in geodesic x = const is 
connected to the comoving one, tc, by the condition ds = {1 — 6p) dtc- 

The key equation (86) is the relativistic Poisson equation and links the 
Laplacian of gravity potential $ with the comoving density perturbation Ssc- 
Equation (88) represents the relativistic Euler equation, or Newton's second 
law. It connects the flow acceleration q with the pressure gradient Spc acting 
on it. 

We have four gauge invariant scalars q , 5 , 6p and S, but the four metric 
potentials A, B, C, and D are not gauge invariant: any two of them can 
arbitrarily be chosen by appropriately selecting the functions X and Y (see 
footnote 21). Thus, in total we end with six independent scalars describing 
density perturbations in the Friedmann model. Gravitational equations (86)- 
(88) impose four constraints on the six potentials of an S mode. Apparently, 
the gravity equations alone are insufficient for describing the dynamics of S 
perturbations. One needs information on the physical state of matter in the 
form of two missing relations (equations of state). However, even without 
imposing constraints on the state of matter, we can derive general evolution 
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equations for the quasi-Friedmann model, similar to the Priedmann equations 
in the homogeneous cosmology. 

In order to describe the geometry of weakly inhomogeneous flows, let us 
use scalar variables, including both orders (the zeroth and first) of pertur- 
bation theory. For the matter, these are the scalars of total density e and 
pressure p, while for the metric they are the scalar scale factor a and the 
Hubble function Hy of the volume expansion. Their relationship is defined 
as 

Hy=l^lu^,,^u^{lnay)^^ (90) 

(the dot above a letter denotes a partial derivative over tc), where we have 
introduced the following factors of medium volume expansion: 

ay = adet^^%6ij-B,ij) = b-{l-q), b = ■ {1 - AB/3) . (91) 

The generalized Priedmann equation (21) for a weakly inhomogeneous Uni- 
verse follows after direct summation of the first equations (82) and (86) and 
grouping terms in the function Hy and s. The other equation (23) follows 
from Eqn (77) with account for the relationship e + ptot = (1 + <5p)(£^ + Pv)- 
Consider the structure of quasi- Hubble flow. According to Eqn (89), the 
proper distance between two neighboring medium elements, separated by 
coordinates Sx\ is given by 

= a {6,, - B,,,) 6x' , det (^^^ = al . (92) 

Differentiating Sri over the proper time s, we obtain the field of paired ve- 
locities of matter motion in the weakly inhomogeneous Universe: 

dSr ■ 1 
5Vi = -^^Hij5r\ Hy^-Hl (93) 

Hj^j = Hq 5ij hij , Hq , hij ^— : (\^ij •) 

cLq a -ij 

where Hij is the matrix of gauge invariant Hubble functions describing the 
recession of matter in space. Only one feature here reminds us of the Hubble 
expansion: the relative recessional velocity for points of the medium is pro- 
portional to the distance between them. Yet these velocities are anisotropic 
and depend on spatial location. 
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B. Dynamics of cosmological scalar perturbations 

If motions of quasi- Hubble flow obey Friedmann equation (21), then the 
cosmological perturbations behave as oscillators. In order to derive the evo- 
lution equation for the adiabatic scalar g, we write the general relationship 
between perturbations of comoving pressure and matter energy density in 
the following form: 

= (5^ + 5p , (94) 

where the function describes the speed squared at which scalar perturba- 
tions propagate in the medium (the speed of sound). For the Pascal media 
(5 = 0), it is as follows: 

P^ = Y.fm Pli , (95) 

m,l 

where is the acoustic matrix of linear medium perturbations!^ The 
scalar of isometric pressure perturbation 5p describes the part of pressure 
which is not related to perturbations of the total energy density 5. Indeed, if 
initially only adiabatic perturbations are present, then 5p = 0, and for 6 = 0, 
one finds 6p = 6p. The relationship between 6p and field variables of medium 
components depends on the equation of state of matter. 
From relationships (88) and (94) it follows that 

6ec ={e + p)6 = a^H (^q - H6p^ , (97) 
2 ^ e + p ^ 7 

Substituting these formulas into the Poisson equation (86), we arrive at 

«V (^q-HSp^ = ja^(3^Aq^. 

^^It connects partial amplitudes of comoving perturbations of pressure and matter den- 
sity, = and the following relationships hold true (see Chapter 6 in book [5] 
for details): 

/™/3.L = M™, /;/3' + /; = E/™/5™^' E/™ = 0- (96) 
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Direct differentiation of this relationship provides the equation for q [2]: 



q+(^H + 2'^q- = l{6„ s) , (98) 

On the left-hand side of equation (98), we have the acoustic d'Alambert op- 
erator for the scalar q describing the general adiabatic density perturbation. 
The right-hand side contains the source of dynamic action of isometric per- 
turbation modes on the evolution of q. For an ideal Pascal media {S = 0), a. 
more compact form of equation follows: 



7a2/3-2 (q' -H5, 



70^ 



Ag = , (99) 



where the prime stands for the derivative over the conformal time r/, and 
H = aH. 

Equations (98) and (99) are valid for a broad variety of media (in par- 
ticular, for fundamental scalar fields). We did not refer to the information 
on the microscopic structure of matter. The only geometric characteris- 
tic of a medium needed for the derivation of key equation (98) is the 
mean velocity squared of propagation of scalar perturbations in the medium 
(95). Neglecting isometric perturbations (for example, in the case of a single 
medium, matter and A-term, and others), 6p = and the equation for q 
acquires a closed form jTUl [H] : 

q+ (3H + 2-]q-(3'^ = 0. (100) 
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C. Eulerian coordinates and the Newton limit 



The Eulerian reference frame (r, y) is uniquely set by the conditions B = 
C = and is of interest, in particular, because the peculiar velocity of matter 

Vpec = = — ^, «^ = (l + ^,i;,i), (101) 

a n. 

where t; = q/if, \I' = $ — 5, is defined relative to this grid. From the 
relationship between Lagrangian and Eulerian coordinates: 

T = tE = t,-v, y = x + S, (102) 

we obtain the metric tensor in the Eulerian representation: 

rfs^ = (1 + 2^)dr2 - a2c^y^ (103) 

where S = — Vi3 is the displacement vector of a medium element relative to 
its initial position. Apparently, the metric in the Eulerian coordinate system 
is independent of the gradients of potential g, and the scalar scale factor 

a = a(r)-(l-<l>) (104) 

fully describes the locally isotropic observer space filled with inhomogeneous- 
matter. Metric (103) is the relativistic limit of the weak field in the Fried- 
mann model, and for the nonrelativistic substance (jiS/^l ^ 1) we obtain 
the Newtonian limit \l/ = $. 

Since the curvature scalar and gravitational potential are small (g ~ $ ^ 
10^^) in the real Universe within the observed structure scale, we can drop 
terms v,'^ , and $ in Eqns (102)- (104), respectively, and introduce the phys- 
ical Eulerian coordinate for the position of medium points r ~ ay. Hence 
we obtain a convenient approximation for the description of the quasilinear 
stage of structure formation in the Universe: 

r = (l + ^)"^(x + S), Vpec = r-i/r = aS. (105) 

Under the additional assumption of the smallness of the speed of sound, from 
equation (100) it follows that g ~ go = Q'o(x). Substituting it into formula 
(105), we obtain the Zel'dovich approximation (see Ref. [Ij for more detail): 

S = -5fVgo, Vpcc = -z/Vgo, 6 = g Ago , $ = 0go, (106) 

where the growth factors are introduced [see Eqn (12)]. They depend only 
on time and the relationship a(j)/g = ?)VL^Hq/2 = const holds true. 
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